ESTCP  199902 


Demo  #2  GPR  UXO  Classification  Results 
for  The  Blossom  Point  Site 

Final  Report 


Matthew  .B.  Higgins\  Chi-Chih  Chenl,  and  Kevin  O’Neill^ 

'Ohio  State  University,  Electroseience  Laboratory 
1320  Kinnear  Rd.,  Columbus,  OH  43212 


2 

US  Army  Corps  of  Engineers 

Engineer  Researeh  and  Development  Center  (ERDC/  CRREL) 
72  Lyme  Rd,  Hanover,  NH  03755 


August  2005 


Distribution  Statement  A;  Approved  for  public  release,  distribution  is  unlimited 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

AUG  2005  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2005  to  00-00-2005 

4.  TITLE  AND  SUBTITLE 

Demo  #2  GPR  UXO  Classification  Results  for  The  Blossom  Point  Site 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Ohio  State  University, Electroscience  Laboratory, 1320  Kinnear 
Rd,Columhus, OH, 43212 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

122 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


Table  of  Contents 


TABLE  OF  CONTENTS . I 

LIST  OF  FIGURES . Ill 

LIST  OF  TABLES . V 

CHAPTER  1  INTRODUCTION . 6 

1 . 1  Radar  System  Description . 2 

1 .2  Measurement  Approach . 3 

CHAPTER  2  IMPROVED  PROCESSING  TECHNIQUES . 6 

2. 1  Improved  Depth  and  Late-Time  Region  Estimation . 6 

2.2  Adaptive  Data  Smoothing  for  Clutter  Suppression . 7 

2.3  Adaptive  Filter  To  Enhance  Target  Resonance . 10 

2.4  Improved  EEF  Estimation  Using  Prony  Coefficients . 10 

2.5  Block  Diagram  of  New  Feature  Extraction  Procedure . 12 

CHAPTER  3  IMPROVED  UXO  CLASSIFICATION  USING  SPATIAL 
FEATURE  AND  SCATTERING  PATTERN . 14 

3.1  Previous  Classification  Problem . 14 

3.2  Examples  from  Known  Targets . 15 

CHAPTER  4  FEATURE  EXTRACTION  AND  UXO  CLASSIFICATION  OF 
BLOSSOM  POINT  TARGETS . 28 

4. 1  Known  Target  Feature  Extraction . 28 

4.2  Offset  Tests  for  Effects  of  Positional  Uncertainty . 32 

4.3  UXO  Classification  Algorithm . 37 

CHAPTER  5  BLIND  TARGETS  CLASSIFICATION  RESULTS  COMPARED 
WITH  GROUND  TRUTH . 40 

5.1  Designation  of  UXO-Like  Items . 40 

5.2  Overall  Classification  Performance . 47 

5.3  Correctly  Classified  UXO-like  Items . 51 

5.4  Missed  UXO-like  Items . 56 

5.5  Clutter  Classified  as  UXO-like  Items . 57 

5.6  Confidence  Level  and  Signal-to-Clutter  Ratio . 70 

CHAPTER  6  ANALYSIS  OF  NOISE  AND  CLUTTER  EFFECTS  ON 
EXTRACTED  PARAMETERS . 71 

6.1  Simulated  Signal  TO  Noise  Ratio  (SNR)  Study  71 

6.2  Measured  Signal  to  Noise  Ratio  (SNR)  80 

6.3  Simulated  Signal  to  Clutter  (SCR)  Study  83 


1 


6.4  Measured  Signal  to  Clutter  Ratio  (SCR)  88 

CHAPTER  7  CONCLUSIONS . 98 

APPENDIX  A  ARCHIVING . 103 

A.l  GPR Data  Files  -  *.CDT . 103 

A.2  Comment  Text  Files  -  *  .txt . 104 

A.3  Processed  Files  -*.mat  Files . 104 

APPENDIX  B  ELECTRICAL  PROPERTIES  OF  THE  BLOSSOM  POINT  SOIL 
. 106 

APPENDIX  C  BLOSSOM  PT  DEMO  LOG  EXCERPT . 110 

BIBLIOGRAPHY . 116 


ii 


List  of  Figures 


Figure  1  Towed  antenna  oriented  parallel/transverse  to  the  motion  direction. 

. 3 

Figure  2  Measurement  passes  performed  at  the  Blossom  Point  Site . 5 

Figure  3  The  background  removed  GPR  data  before  applying  adaptive  data 

SMOOTHING . 7 

Figure  4  The  background  removed  GPR  data  after  applying  adaptive  data 

SMOOTHING . 8 

Figure  5  Comparison  of  GPR  data  before  (a)  and  after  (b)  adaptive  data 

SMOOTHING.  The  target  is  a  vertical  UXO  (Blossom  Point  Known  Target  C).  9 
Figure  6  Comparison  of  different  method  of  extracting  ELF  feature  in  the 

PRESENCE  OF  DIFFERENT  NOISE  LEVEL .  11 

Figure  7  Block  diagram  of  the  improved  UXO  feature  extraction  algorithm.  1 3 
Figure  8  Different  scenarios  of  linear  and  non-linear  objects  beneath  the  HFB 

ANTENNA .  15 

Figure  9  Pictures  of  known  targets  at  Blossom  Point  site . 16 

Figure  10  Preprocessed  GPR  data  for  Target  B  (horizontal  solid  steel 

CYLINDER,  1.5"  X  6") .  18 

Figure  1 1  Extracted  ELE  and  ETO  for  Target  B  (horizontal  solid  steel 

CYLINDER,  1 .5"  X  6")  AS  A  FUNCTION  OF  ANTENNA  POSITION .  19 

Eigure  12.  Preprocessed  GPR  data  for  Target  M  (vertical  MK23) . 20 

Eigure  1 3 .  Extracted  ELE  and  ETO  for  Target  M  (vertical  MK23)  as  a  function 

OF  ANTENNA  POSITION . 21 

Eigure  14  Preprocessed  GPR  data  for  Target  G  (16lb  shotput) . 22 

Eigure  1 5 .  Extracted  ELE  and  ETO  for  Target  G  ( 1  6lb  shotput)  as  a  function 

OF  ANTENNA  POSITION . 23 

Eigure  16  Preprocessed  GPR  data  for  Target  H  (horizontal  plate, 

4"  X  4"  X  0.25") . 24 

Eigure  1 7  Extracted  ELE  and  ETO  for  Target  H  (horizontal  plate,  4"  x  4"  x 

0.25")  AS  A  FUNCTION  OF  ANTENNA  POSITION . 25 

Eigure  18.  Preprocessed  GPR  data  for  Target  J  (horizontal  plate,  8"  x  8"  x 

0.25") . 26 

Eigure  19.  Extracted  ELE  and  ETO  for  Target  J  (horizontal  plate,  8"  x  8"  x 

0.25")  AS  A  FUNCTION  OF  ANTENNA  POSITION . 27 

Eigure  20.  Energy  plot  of  the  measured  S  1 1  data  for  the  vertical  MK23 

(Known  Target  M) . 29 

Eigure  2 1  Bipolar  plot  of  the  measured  S  1 1  data  for  the  vertical  MK23  (Known 

Target  M) . 30 

Eigure  22.  Comparisons  of  true  and  estimated  known  targets  lengths  (top) 

AND  DEPTHS  (BOTTOM),  SORTED  BY  MAGNITUDE  ORDER . 31 

Eigure  23 .  A  north-directed  scan  with  its  scan  center  offset  from  target 

CENTER  BY  1 0  INCHES  FROM  BOTH  EAST  AND  NORTH  DIRECTION . 33 


iii 


Figure  24  Classification  features  and  position-time  plot  for  target  B  1 5 , 

OBTAINED  FORM  THE  10-INCH  (14  INCH  DIAGONAL)  OFFSET  SURVEY  PATH . 33 

Figure  25  An  East-West  scan  with  its  scan  center  offset  from  target  center 

BY  20  INCHES  FROM  BOTH  EAST  AND  NORTH  DIRECTION . 34 

Figure  26  Classification  features  and  position-time  signal  plot,  for  target 

B15,  OBTAINED  FROM  THE  20-INCH  OFFSET  SURVEY  PATH . 34 

Figure  27  Classification  features  for  a  vertical  MK23  bomb  obtained  form 

THE  10-INCH  OFFSET  CONFIGURATION . 35 

Figure  28  Classification  features  and  position-time  signal  plot  for  a  vertical 
MK23  BOMB  OBTAINED  FROM  A  20-INCH  OFFSET  (28  INCH  DIAGONAL)  SURVEY  PATH. 

. 36 

Figure  29.  Improved  UXO  classification  flow  chart . 39 

Figure  30.  Pictures  of  the  old  Blossom  Point  items  in  lanes  A~E . 42 

Figure  3 1  Pictures  of  the  new  Blossom  Point  items  in  lanes  F~1 . 46 

Figure  32.  The  classification  performance  for  all  targets  buried  at  the 

Blossom  Point  site  considered  together  (old  ones  lanes  A~E  and  new  ones 

LANES  E~l),  INCLUDING  CALIBRATION  TARGETS  TREATED  AS  UNKNOWNS . 49 

Eigure  33  Classification  performance  for  the  old  targets  alone  (A~E  lanes) 

BURIED  AT  THE  BLOSSOM  POINT  SITE,  INCLUDING  CALIBRATION  TARGETS  TREATED  AS 

UNKNOWNS . 49 

Eigure  34.  The  classification  performance  for  the  new  targets  (E~1  lanes) 

BURIED  AT  THE  BLOSSOM  POINT  SITE . 50 

Eigure  35 .  The  classification  performance  for  all  lanes  excluding  the 

CALIBRATION  TARGETS . 50 

Eigure  36.  The  classification  performance  for  the  old  lanes  excluding  the 

CALIBRATION  TARGETS . 51 

Eigure  37.  Comparison  of  estimated  depth  and  true  depth  of  correctly 

CLASSIFIED  UXO-LIKE  ITEMS . 54 

Eigure  3  8 .  Comparison  of  estimated  length  and  true  length  of  correctly 

CLASSIFIED  UXO-LIKE  ITEMS . 55 

Eigure  39.  Comparison  of  estimated  orientation  and  true  orientation  of 

CORRECTLY  CLASSIFIED  UXO-LIKE  ITEMS . 55 

Eigure  40.  Modified  classification  rules  to  include  features  from  vertical 

PLATES . 59 

Eigure  4 1 .  Modified  classification  rules  to  include  features  from  broken 

LOOPS  AND  VERTICAL  PLATES . 60 

Eigure  42.  Classification  features  for  Item  E8  (fence  post  cap)  from  survey 

LINE  PARALLEL  TO  TARGET  ORIENTATION,  FOR  FENCE  POST  CAP . 61 

Eigure  43  Classification  features  for  Item  A9  (vertical  plate)  obtained  from  a 

PASS  ORTHOGONAL  TO  TOP  EDGE  DIRECTION  OF  THE  PLATE . 62 

Eigure  44  Classification  features  for  Item  H  1 1  (crumpled  box  fin)  from  an 

OFFSET  CUT . 63 

Eigure  45  Classification  features  for  Item  C9  (horse  shoe) . 64 

Eigure  46  Classification  features  for  Item  D  1  (M3  8  box  fin  piece)  from  pass 

PARALLEL  TO  TARGET  ORIENTATION . 65 

Eigure  47  Classification  features  for  Item  D8  (M3 8  box  fin) . 66 


IV 


Figure  48  Classification  features  for  Item  H8  (horizontal  8”  dia.  steel  disk) 

. 67 

Figure  49.  Classification  features  for  Item  H  1 2  (box  fin  fragment)  from  pass 

PARALLEL  TO  TARGET  ORIENTATION . 68 

Figure  50  Classification  features  for  Item  B  1 3  (vertical  3”x6”  cyl.)  from 

SURVEY  PASS  PARALLEL  TO  TARGET  ORIENTATION . 69 

Figure  5 1  Classification  features  for  Item  C  1 0  (vertical  4”x4”  plate)  from 

PASS  PERPENDICULAR  TO  TARGET  ORIENTATION . 69 

Figure  52  Figure  53 . 70 

Figure  54  Conceptual  smart  antenna  design  with  narrow  beam  patterns  (red) 

AND  CONTROLLABLE  BEAM  DIRECTION,  SHOWN  ORIENTED  IN  SELECTED  DIRECTIONS. 


Figure  55  GPR  data  collected  parallel  to  the  UXO  (2.75”,  Position  19) 

ORIENTATION  THAT  HAS  A  30-DEGREE  DEPRESSION  ANGLE .  lOI 

Figure  56  GPR  data  collected  parallel  to  the  UXO  (8 1  mm.  Position  1 1 3) 

ORIENTATION  THAT  HAS  A  30-DEGREE  DEPRESSION  ANGLE .  102 

Figure  5  7  Stronger  scattering  for  a  tilted  UXO  from  offset  position . 102 

Figure  58  OSU/ESL  soil  probe  for  permittivity  and  conductivity  measurement 

AT  40  MHz  AND  60  MHz . 106 

Figure  59  Soil  electrical  property  measured  on  12/10/2000  at  40  MHz . 107 

Figure  60  Soil  electrical  property  measured  on  12/10/2000  at  60  MHz . 107 

Figure  61  Soil  electrical  property  measured  on  12/1 1/2000  at  40  MHz . 108 

Figure  62  Soil  electrical  property  measured  on  12/1 1/2000  at  60  MHz . 108 

Figure  63  Soil  electrical  property  measured  on  12/12/2000  at  40  MHz . 109 

Figure  64  Soil  electrical  property  measured  on  12/12/2000  at  60  MHz . 109 

Figure  65  Location  designation  of  the  Blossom  Point  site.  Only  lanes  AM  were 

MEASURED  BY  OSU/ESL  GPR .  1 12 

Figure  66  Target  ID  of  the  new  GPR  targets  buried  in  lanes  FM . 114 

List  of  Tables 

Table  1  Descriptions  of  Known  Targets . 16 

Table  2  Locations,  Depths  and  Orientations  of  Know  Targets . 16 

Table  3  Features  Extracted  for  Known  Targets . 30 

Table  4  UXO-Like/Clutter  Designation  of  Old  Blossom  Point  Items  (Lanes  A~E) . 41 

Table  5  UXO-Eike/Clutter  Designation  of  New  Blossom  Point  Items  (Eanes  F~I) . 43 

Table  6  Eist  of  correctly  classified  UXO-like  items . 53 

Table  7  Eist  of  missed  UXO-like  items . 56 

Table  8  Eist  of  clutter  classified  as  UXO-like  items . 57 

Table  9  Format  for  multiple-position  GPR  data . 103 

Table  1 0  Data  format  for  a  single  frequency  at  a  single  position . 104 

Table  1 1  The  target  ID  (Ground  Truth)  of  the  targets  buried  in  Eanes  A~E . 113 

Table  12  Description  of  the  targets  buried  in  lanes  A~E . 1 14 

Table  13  Description  of  target  buried  in  lanes  F~I . 115 


V 


Chapter  1  Introduction 


An  ultra-wide  bandwidth,  full-polarimetric  ground  penetrating  radar  (GPR)  system  operating 
from  10  MHz  to  810  MHz  has  been  applied  to  the  classifieation  of  unexploded  ordnance  (UXO). 
Late -time  radar  signatures  such  as  natural  resonance  and  polarization  were  utilized  to  separate 
false  alarm  objects  from  UXO-like  items  with  elongated  bodies  with  length- to-diameter  ratios 
greater  than  three.  Additional  features  such  as  length,  depth  and  orientation  were  also  extracted 
from  the  data.  Good  classification  performance  was  obtained  from  the  previous  Tyndall  Test  [1] 
for  UXO-like  items  with  inclination  angles  less  than  45  degrees  [2],  where  inclination  is 
measured  relative  to  the  approximately  horizontal  ground  surface.  That  particular  test  utilized 
only  features  extracted  from  a  single  radar  position,  approximately  directly  above  the  target.  It 
was  later  found  that  this  single-position  approach  contributed  to  many  false  alarms  and  the 
incapability  of  classifying  UXO’s  with  large  inclination  angles,  such  as  vertical  UXO’s. 
Subsequently,  new  measurement  approaches  and  additional  radar  features  were  developed  to 
reduce  the  false  alarm  rate  and  missed-UXO’s. 

A  second  field  test  (“Demo  #2),  conducted  at  the  Blossom  Point  UXO  Site  [3,4],  is  the 
subject  of  this  report.  An  improved  system  configuration  and  measurement  approach  was  used 
to  collect  GPR  data  at  multiple  positions  by  moving  the  antenna  along  a  straight  line  centered  at 
the  target  spot.  Other  than  14  known  targets,  there  were  78  blind  target  positions  marked,  where 
either  UXO’s  or  clutter  items  had  been  emplaced  without  the  test  team  knowing  the  target 
identities.  Multiple  passes  were  performed  along  different  orientations  through  each  target 
position.  The  new  approach  provides  spatial  distributions  of  scattering  response,  including  both 
early- time  behavior  as  well  as  late-time  UXO  signatures.  The  combination  of  these  features 
improves  the  classification  of  UXO’s  with  large  inclination  angles  due  to  the  additional  oblique 
observation  angle  with  respect  to  the  target.  False  alarm  items  caused  by  offset  shallow  objects 
can  now  be  discriminated  using  the  additional  time-position  scattering  characteristics.  Overall, 
the  current  improved  approach  should  increase  the  success  rate  of  UXO  classification  and 
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decrease  the  false  alarm  rate.  The  new  approach  also  eliminates  the  problem  of  position 
uncertainty  along  the  scan  direction  because  it  inherently  samples  over  a  structured  variety  of 
positions  and  measurement  pass  orientations. 

This  report  summarizes  the  new  processing  and  classification  algorithm,  and  presents  the 
classification  results  for  all  the  targets  measured  at  the  Blossom  Point  site. 


1.1  Radar  System  Description 

The  same  radar  system  used  in  the  Tyndall  Demonstration  was  used  again  at  the  Blossom 
Point  site  except  that  the  new  frequency  range  was  from  10  MHz  to  410  MHz  at  2  MHz 
increments.  Notice  that  the  upper  frequency  and  the  increment  have  been  changed  from  810 
MHz  and  1  MHz,  respectively,  used  in  the  earlier  Demo.  This  is  because  the  Blossom  Point  soil 
was  quite  wet,  with  high  dielectric  constant.  Such  conditions  render  the  higher  frequencies 
useless  because  they  are  absorbed  by  the  moist  soil;  further,  under  these  conditions  the  resonance 
phenomena  of  greatest  interest  appear  in  a  lower  frequency  range.  The  larger  frequency 
increments  were  used  to  speed  up  data  collection.  The  radar  includes  a  commercial  RF  vector 
network  analyzer  (HP8712ET)  using  a  dielectric  loaded  horn-fed  bow-tie  (HFB)  antenna.  It  is 
set  up  for  collecting  both  co-polarized  and  cross  polarized  returns.  Relative  to  the  previously 
used  configuration,  the  antenna  orientation  was  rotated  45  degrees  such  that  the  two  dipole- 
elements  are  parallel  and  transverse  to  the  direction  of  scan  as  shown  in  Figure  1 .  The  co¬ 
polarized  returns  are  Sn  (transmission  and  reflection  polarization  along  the  vehicle  path),  S22 
(polarization  perpendicular  to  the  vehicle  path),  and  S12  and  S21  are  the  two  cross-polarizations. 
Assuming  that  the  target  positioning  information  is  approximately  correct,  the  new  antenna 
orientation  ensures  the  target  remains  in  the  principal  plane,  which  has  purer  radiation 
polarization  during  the  position  scan.  The  antenna  was  towed  behind  an  ATV  such  that  it  could 
be  lifted  off  the  ground  when  necessary. 

The  radar  control  and  data  collection  was  done  from  a  laptop  computer  via  special  purpose 
software  developed  during  this  project.  Both  frequency-domain  and  time-domain  data  are 
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displayed  in  real  time  so  that  the  operator  can  monitor  the  status  of  the  operation.  Although  on¬ 
site  processing  can  be  performed  immediately  after  each  pass,  the  actual  processing  was 
performed  overnight  to  maximize  the  data  collection  rate. 


Figure  1  Towed  antenna  oriented  parallel/transverse  to  the  motion  direction, 

1.2  Measurement  Approach 

The  focus  of  this  project  is  on  discrimination,  not  detection.  That  is,  we  assume  that  an 
earlier  pass  with  a  magnetometer  or  metal  detector  has  produced  position  information  for 
anomalies  which  we  will  investigate  further  to  produce  a  UXO-like  vs  non-UXO  classification. 
Our  initial  pass  through  each  designated  target  spot  was  conducted  at  an  arbitrarily  chosen 
direction  unless  available  magnetic  sensor  data  provided  a  strong  suggestion  of  the  orientation  of 
the  target.  Each  pass  length  was  about  10  feet  long  with  a  sample  spacing  of  3  inches.  At  each 
sampling  location,  full-polarimetric  data  (both  amplitude  and  phase,  in  two  co-polar  and  one 
cross  polar  orientation)  was  collected  from  the  two  orthogonal  antenna  elements  from  10  MHz  to 
410  MHz  at  a  2  MHz  increment.  This  information  can  be  manipulated  mathematically  to 
produce  equivalent  co-  and  cross-polarized  returns  that  would  be  expected  from  any  other 
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polarization  orientation.  This  full-polarimetric  ultra-wideband  measurement,  not  available  in  any 
other  GPR  systems  we  know  of,  offers  a  unique  eapability  for  deriving  signal  parameters  whieh 
allow  the  inferenee  that  an  elongated  target  is  present.  Further,  the  length  of  the  target  may  be 
estimated  as  well,  based  on  a  eoneomitant  estimation  of  soil  dieleetrie  properties. 

The  data  from  a  days  surveying  is  easily  stored  on  deviees  attaehed  to  the  rig;  it  is  then 
proeessed  overnight  rather  than  in  real  time,  in  order  to  speed  up  the  data  eolleetion  rate.  After 
an  initial  pass  for  eaeh  target,  a  single  additional  pass  orthogonal  to  the  initial  pass  was  taken  for 
any  target  that  showed  poor  resonanee  and  preferred  polarimetrie  orientation.  A  preferred 
polarization  that  produees  late-time  resonanee  suggests  the  orientation  of  an  elongated  target. 
Absenee  of  resonanee  offers  no  preferred  orientation  to  guide  subsequent  passes,  and  further 
suggests  that  the  objeet  is  not  UXO-like.  For  targets  that  showed  signifieant  resonanee  and 
linearity  features  from  the  overnight  proeessing,  two  additional  passes  were  made  in  the 
direetions  parallel  and  transverse  to  the  estimated  target  orientation  extraeted  from  the  initial- 
pass  data.  The  idea  was  of  these  latter  two  pass  orientations  was  to  get  views  of  the  target  with 
polarization  optimally  aligned  with  its  prineipal  direetions,  rather  than  having  to  infer  data  in 
those  direetions.  Figure  2  indieates  all  the  passes  performed  at  the  Blossom  Point  site.  To  test 
for  the  effeets  of  positional  uneertainty,  additional  passes  were  performed  for  several  ealibration 
targets  (Rows  IS-IS)  with  the  introduetion  of  lO-ineh  and  20-inoh  position  offsets  along  the 
225 -degree  direetion. 
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Direction  of  Scans  for  Blossom  Point  Measurements 


ABCDEFGHI 


Figure  2  Measurement  passes  performed  at  the  Blossom  Point  Site, 
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Chapter  2  Improved  Processing  Techniques 


This  chapter  discusses  several  new  data  processing  techniques  developed  to  significantly 
improve  the  accuracy  of  the  extracted  features. 


2.1  Improved  Depth  and  Late-Time  Region  Estimation 

The  full-polarimetric,  broadband  GPR  data  collected  from  different  positions  and  passes  for 
each  target  was  processed  in  a  similar  technique  used  in  the  Demo  I  for  each  position.  A 
background  response  for  each  channel  for  each  pass  was  generated  from  the  ensemble  average  of 
all  positions  in  that  pass.  For  background  removal  (e.g.  layering)  this  data  was  then  subtracted 
from  data  at  each  antenna  position,  before  the  feature  extraction  procedure. 

During  the  Demo  I,  the  target  portion  of  response  was  based  on  a  single  time  waveform 
collected  at  a  single,  hopefully  optimal  position.  The  optimal  position  was  often  not  so  clear. 

This  can  occur  when  the  necessary  signal  level  is  not  much  greater  than  the  clutter  level,  as  in  the 
cases  of  deep  target  where  the  target  response  is  weak,  or  for  a  shallow  small  target  where 
surface  clutter  is  significant.  An  improved  method  of  selecting  the  onset  of  the  late-time 
response  associated  with  the  target  is  made  possible  by  the  multiple  position  data.  With 
additional  position  scans,  the  target  produces  a  hyperbolic  arc,  as  appears  in  most  conventional 
GPR  data  plots  showing  signal  magnitude  through  time  (vertical  axis)  vs  antenna  position 
(horizontal  axis).  We  designate  this  as  a  time-position  plot.  The  inhomogeneous  layers  do  not 
have  this  temporal-spatial  feature  as  demonstrated  in  Figure  3  for  a  UXO  buried  at  75  cm  depth. 
Observing  the  signal  travel  time  to  the  peak  of  the  arc  also  improves  the  accuracy  of  depth 
estimation,  relative  to  using  data  from  a  single  position. 


6 


Arcs  such  as  that  in  the  figure  appear  because,  as  the  antenna  moves  further  away  from  the 
position  directly  above  the  target,  the  signal  travel  time  to  and  from  the  target  becomes  longer. 
Thus  for  offset  antenna  positions  the  target  response  appears  later  in  time.  The  relatively  clear 
are  in  the  figure  is  an  early-time  response,  resulting  from  the  first  refieetion  of  the  signal  from 
the  target.  While  they  are  not  easily  visible  without  further  proeessing  (see  next  section),  this 
initial  arc  is  followed  by  a  succession  of  later,  weaker  ares,  representing  reverberations  of  a  sort, 
i.e.  late  time  resonanee.  Using  a  graphieally  interactive  routine,  the  user  of  the  data  analysis 
program  simply  designates  a  window  on  the  time  position  plot,  delimiting  a  region  clearly  below 
the  first  response  represented  by  the  hyperbolic  arc.  The  exact  boundaries  of  this  window  are  not 
crueial,  beeause  more  detailed  proeessing  of  the  signals  within  this  window  indieates  where 
resonance  oeeurs. 


Blossom  Point  Known  Target  F  (0.75m) 


-60  -40  -20  0  20  40  60 

Antenna  Position  (inches) 

Figure  3  Time-  position  plot  of  GPR  data,  with  background  removed,  before  application 
of  adaptive  data  smoothing. 


2.2  Adaptive  Data  Smoothing  for  Clutter  Suppression 
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A  new  special  algorithm  that  performs  data  smoothing  along  these  arcs  has  been 
implemented  to  greatly  improve  the  signal-to-clutter  ratio.  In  addition  to  the  arcs  associated  with 
the  desired  target,  there  are  also  other  arcs  with  different  intensities  and  curvatures  resulting  from 
surface  and  subsurface  clutter.  This  phenomenon  can  be  clearly  observed  in  Figure  3.  These 
undesirable  arcs  interfere  with  the  target  responses  and  cause  errors  and  instability  in  feature 
extraction.  By  applying  a  smoothing  window  that  runs  along  the  target-arc  curves  in  the  whole 
data  region,  such  interference  can  be  reduced.  This  is  demonstrated  by  comparing  the  smoothed 
data.  Figure  4,  with  the  data  before  smoothing.  Figure  3.  The  arc  shape  for  processing  is 
generated  based  on  a  few  points  selected  by  the  operator.  Again,  precise  selection  of  the  points  is 
not  crucial.  Figure  4  shows  that,  after  the  smoothing,  clutter  has  been  reduced  effectively  and  late 
time  arcs  are  now  visible,  approximately  parallel  to  the  original  early  time  signal.  Another 
example  is  given  for  a  vertically  buried  UXO  (Blossom  Target  C,  Figure  5),  which 
characteristically  produces  a  weaker  peak  and  stronger  sides  in  the  arc.  Note  the  improved  data 
quality  in  the  late -time  region  below  the  initial  arc.  A  group  of  arcs  associated  with  the  trench  is 
also  observed  near  the  surface  (-7  ns  position). 
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Figure  4  The  background-removed  GPR  data  after  application  of  adaptive  data  smoothing 


Blossom  Point  Known  Target  F  (0.75m) 


-50  0  50 

Antenna  Position  (inches) 
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0 


Blossom  Point  Known  Target  C 


Figure  5  Comparison  of  GPR  data  before  (a)  and  after  (b)  adaptive  data  smoothing.  The 
target  is  a  vertical  UXO  (Blossom  Point  Known  Target  C), 
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2.3  Adaptive  Filter  To  Enhance  Target  Resonance 


An  adaptive  bandpass  filter  using  the  Kaiser  window  funetion  is  applied  to  the  late -time 
speetrum  to  seleet  the  resonant  peak.  This  filter  is  adaptive  beeause  its  eenter  frequeney  and 
bandwidth  are  determined  aeeording  to  the  target’s  resonant  speetrum.  The  Kaiser  window  was 
ehosen  beeause  it  produees  low  sidelobe  levels  in  the  time  domain,  when  the  frequency  domain 
data  is  inverted.  Using  a  convenient  graphical  software  interface  to  select  the  dominant  spectral 
peak  and  its  approximate  width,  the  operator  determines  the  filter’s  center  frequency  and  relative 
bandwidth.  This  is  normally  obtained  from  the  late-time  spectrum  from  a  channel  that  produced 
the  strongest  target  responses. 


2.4  Improved  ELF  Estimation  Using  Prony 
Coefficients 

Our  recent  studies  show  that  improved  accuracy  of  the  ELF  estimation  under  the  influence  of 
random  noise  could  be  obtained  by  utilizing  the  autocorrelation  process  and  Prony  model 
coefficients.  Previous  processing  extracts  the  ELF  and  other  linear  feature  parameters  from  each 
temporal  position  in  the  selected  late-time  response.  Then  a  single  parameter  was  obtained  by 
averaging  over  the  selected  late-time  region  to  suppress  the  noise  effect.  First,  the  new 
processing  calculates  the  autocorrelation  responses  of  the  selected  late -time  responses.  The 
Prony  model  is  applied  to  fit  the  autocorrelation  response  of  each  channel  (***at  some  point 
identify  these  with  “channels”  Sn,  S21  and  S22)  with  complex  exponentials.  The  magnitude  of  the 
late-time  response  in  each  channel  is  now  represented  by  a  single  dominant  Prony  coefficient. 
The  two  eigenvalues  and  eigenfunctions  are  then  calculated  from  the  Matrix  form  by  the  Prony 
coefficients  from  the  three  channels.  Recall  that  these  eigenvalues  and  eigenfunctions  are  used 
for  ELF  and  ETO  estimation.  Notice  that  the  eigenvalues  and  eigenfunctions  only  have  to  be 
solved  once  for  one  waveform.  The  previous  approach  required  them  to  be  solved  from  every 
time  position  in  the  late-time  region.  The  final  two  eigenvalues  and  eigenfunctions  were  then 
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obtained  from  the  temporally  averaged  values.  The  magnitudes  of  the  dominant  resonanee 
eomponents  in  each  channel  are  then  used  to  form  a  2x2  scattering  matrix  from  which  ELF,  ETO 
and  DEN  are  estimated.  It  should  be  noted  that,  in  the  previous  processing,  Prony  processing 
was  applied  only  to  the  strongest  response  channel  for  extracting  the  resonance  information. 

This  new  approach  is  more  effective  than  the  previous  averaging  approach  in  suppressing 
random  noise.  Figure  6  compares  the  extracted  ELF  in  the  presence  of  different  noise  levels  from 
the  previous  approach  to  that  obtained  from  the  new.  The  noise  level  is  expressed  as  SNR,  i.e. 
signal  to  noise  ratio,  usually  measured  in  dB.  More  precise  definitions  and  analyses  of  SNR  and 
signal  to  clutter  ratio  (SCR)  are  to  found  below  in  Chapter  6,  where  effects  of  those  two  types  of 
signal  contamination  on  extracted  parameters  are  pursued.  Here  the  signal  is  assumed  to  be  a 
linearly  polarized  damped  resonance  characterized  by  cos(0.3t)exp'°  ‘’^'  with  Gaussian  distributed 
random  noise  added.  These  results  were  obtained  using  500  noise  realizations.  The  improvement 
of  the  new  Prony  approach  is  quite  significant  for  SNR  greater  than  10  dB.  For  SNR  less  than  10 
dB,  almost  all  processing  methods  will  have  difficulty  in  obtaining  the  correct  ELF,  i.e.  1. 


500  Simutations 


Figure  6  Comparison  of  different  method  of  extracting  ELF  feature  in  the  presence  of 
different  noise  level. 
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2.5  Block  Diagram  of  New  Feature  Extraction 
Procedure 


Figure  7  shows  the  block  diagram  of  the  improved  processing  procedure.  The  details  of  the 
feature  extraction  block  are  similar  to  those  used  in  Demo  I  and  can  be  found  in  [1,2],  The  final 
output  contains  the  extracted  ELF,  DEN,  ETO  and  CNR  features  as  a  function  of  antenna 
position.  Several  new  processing  blocks  (shadowed)  have  been  discussed  in  the  previous 
paragraph  and  are  introduced  to  improve  the  accuracy  in  the  presence  of  noise  and  clutter. 
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Figure  7  Block  diagram  of  the  improved  UXO  feature  extraction  algorithm. 
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Chapter  3  Improved  UXO  Classification  Using 

Spatial  Feature  and  Scattering  Pattern 

3.1  Previous  Classification  Problem 

Many  false  alarms  were  caused  by  unbalanced  excitation  due  to  the  target  offset  near  the 
fully-polarimetric  antenna  as  shown  in  Figure  8(a)  and  (b).  This  causes  the  extracted 
polarization  feature  to  resemble  that  obtained  from  a  truly  linear  object,  i.e.  an  object  like  a 
UXO,  as  shown  in  Figure  8(c).  Utilizing  the  spatial  distribution  of  features  obtained  from 
multiple-position  data,  the  non-linear  object  shows  low  linearity  as  the  antenna  is  moved  close  to 
the  top  of  the  target.  This  makes  it  easy  to  be  separated  from  a  UXO-like  target  with  a  small 
inclination  angle.  However,  from  this  feature  alone,  one  still  could  not  discriminate  the  non¬ 
linear  object  from  a  steeply  inclined  UXO  as  illustrated  in  Figure  8(d).  A  similar  null  of  linearity 
will  be  observed  from  this  type  of  UXO  although  the  good  linearity  observed  at  offset  positions 
may  result  from  oblique  excitations.  Additional  early-time  scattering  features  due  to  the  two 
scattering  mechanisms  associated  with  the  nose  and  tail  of  the  UXO  can  sometimes  be  used  to 
separate  it  from  a  small  non-linear  object  that  contains  only  one  scattering  point.  Sometimes  a 
large  horizontal  plate  also  has  multiple  scattering  edges,  its  strong  reflection  from  the  flat  surface 
makes  it  easy  to  distinguish  from  a  vertical  UXO.  These  will  be  demonstrated  by  actual 
measurement  examples. 
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(a)  top  view  w/  a  non-linear  target  beneath  the  ^  non-linear  target  beneath  the  antenna 

antenna 


(e)  a  horizontal  linear  target  beneath  the  antenna  (d)  a  vertieal  linear  target  beneath  the  antenna 

Figure  8  Different  scenarios  of  linear  and  non-linear  objects  beneath  the  HFB  antenna. 


3.2  Examples  from  Known  Targets 

A  set  of  known  targets  was  available  for  tuning  the  system  and  processing.  The  pictures  and 
descriptions  of  these  targets  are  shown  in  Figure  9  and  Table  1,  respectively.  Table  2  also  shows 
their  locations,  depths  and  orientations. 
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Figure  9  Pictures  of  known  targets  at  Blossom  Point  site. 


Table  1  Descriptions  of  Known  Targets 


Item  ID 

Description 

Item  ID 

Description 

A 

1 14"  X  3"  solid  steel 
cylinder 

G 

16  lb  steel  shotput 

B 

1  V-i  X  6"  solid  steel 
cylinder 

H 

4"  X  4"  X  yv  steel  plate 

C 

r/2"  X  12"  solid  steel 
cylinder 

J 

8"  X  8"  X  yv  steel  plate 

D 

3"  X  6"  solid  steel  cylinder 

K 

iy2"x6"xy4"  steel  plate 

E 

3"x  12"  solid  steel 
cylinder 

L 

3"  X  12"  X  yv  steel  plate 

F 

3"  X  24"  solid  steel 
cylinder 

M 

Table  2  Locations,  Depths  and  Orientations  of  Know  Targets 


Row/Column 

A 

B 

C 

D 

E 

15 

16-lb  shotput 
25  cm 

Item  L 

25  cm/  0°  /  0° 

Item  B 

20  cm/  0°  /  0  ° 

Item  A 

5  cm/  0°  /  0° 

Item  K 

8  cm/  0°  /  90° 

14 

Item  H 

5  cm/  45°  /  0° 

Item  H 

25  cm/  0°  / 
90° 

Item  J 

25  cm/  45°  / 

0° 

Item  C 

35  cm/  0°  / 
90° 

Item  F 

75  cm/  0°  /  0° 
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Item  P 

75  cm/  0°  / 
90° 

Item  D 

10  cm/  0°/ 
90° 

Item  E 

50  cm/  0°  / 
90° 

Item  N 

35  cm/  0°  / 
90° 

empty 
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The  relative  permittivity  and  eonductivity  of  the  soil  at  60  MHz  are  shown  in  Appendix  B  as  a 
funetion  of  depth.  Figure  10  plot  the  delay  time  vs.  antenna  position  plot  of  the  energy  of  the 
measured  fully-polarimetric  data  for  a  horizontal  UXO-like  Object  (Target  B).  The  grayscale 
indicates  the  amplitude  of  the  scattered  fields.  Notice  that  S22  channel  contains  most  significant 
response  because  antenna  2  is  parallel  to  the  UXO  axis.  Figure  1 1  shows  the  estimated  linear 
factor  (ELF)  and  orientation  (ETO)  as  a  function  of  position.  The  large  ELE  values  when  the 
antenna  is  near  the  target  location  clearly  indicated  a  linear  object,  i.e.  an  ElXO-like  target. 
Accurate  orientation  estimations  are  also  obtained. 

Eigure  12  and  Eigure  13  plot  the  measured  and  processed  results  for  a  buried  vertical  MK23 
(Target  M).  Notice  the  low  ELE  region  near  the  center  because  the  EIXO  has  a  rotational 
symmetry  when  observed  directly  above.  Sn  shows  the  dominant  responses  and  shows  maximal 
response  when  the  antenna  is  offset  from  the  target.  This  is  due  to  the  stronger  excitation  of 
resonant  modes  of  the  EIXO  when  incidence  is  from  the  side.  In  general,  we  will  see 
insignificant  scattering  for  an  EJXO  when  the  direction  of  incidence  is  directly  along  the  EIXO 
axis,  within  our  operational  frequency  range,  except  for  an  extremely  large  EJXO.  Two  scattering 
arcs  that  have  different  curvatures  are  clearly  observed  between  position  15  and  40  in  Eigure 
12(a).  These  are  generated  from  the  nose  and  tail  of  the  EJXO.  The  smaller  curvature  one  is 
generated  form  the  bottom  of  the  EJXO  and  the  larger  curvature  one  is  generated  by  the  top  of  the 
EJXO.  The  estimation  orientations  are  aligned  with  the  orientation  of  the  antenna  1  . 

Eigure  14  and  Eigure  15  plot  the  measured  and  processed  results  for  a  shallow  conducting 
sphere  (Target  G).  The  single  arc  pattern  indicates  only  one  scattering  point.  This  is  different 
from  what  observed  in  the  previous  vertical  EJXO  case.  Notice  that  the  ELE  peaks  up  in  two 
regions  offset  from  the  center.  These  regions  correspond  to  the  portions  of  the  antenna  1  that  are 
closest  to  the  ground  surface  as  illustrated  in  Eigure  8(a).  It  is  not  surprising  to  find  that  the 
estimation  orientations  are  aligned  with  the  orientation  of  the  antenna  1  (90  degrees  in  this  case). 
Similar  results  can  be  obtained  from  a  small  plate  (Target  H)  as  shown  in  Eigure  16  and  Eigure 
17. 

Eigure  18  and  Eigure  19  plot  the  measured  and  processed  results  for  a  large,  horizontal 
conducting  plate  (Target  J).  Similar  to  the  sphere  case,  the  ELE  peaks  up  in  two  regions  offset 
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from  the  center.  Notice  the  strong  reflection  observed  in  middle  of  both  Sn  and  S22  channels. 
This  is  distinctive  from  a  sphere  or  a  vertical  UXO  except  for  a  very  large  nose-down  UXO  that 
has  a  flat  end. 


•50 


■60 -  - iQ 

-60  -40  -20  0  20  40  60 

Position  (inches) 

(a)  Sii  data  for  Target  B 


S21  Hem  ^  .jq-3 


Position  (inches) 


(b)  S21  data  for  Target  B 


S22  Item  g'- ^  ^q-3 


(e)  S22  data  for  Target  B 


Figure  10  Preprocessed  GPR  data  for  Target  B  (horizontal  solid  steel  cylinder,  1,5"  x  6"), 
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Frequency  Domain  Estimated  Linear  Factor 


1 


fELF  vs  Position  for  C15  354deg 


(a)  ELF  data  for  Target  B 


ETO  vs  Position  for  C15  354de9 


(b)  ETO  for  Target  H 

Figure  11  Extracted  ELF  and  ETO  for  Target  B  (horizontal  solid  steel  cylinder,  1,5"  x 
6")  as  a  function  of  antenna  position. 
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S11  ItemM 


X  10 


S21  Item  M*- 


-60  -40 


20  40 

Position  (Inches) 


-60  -40 


(a)  Si  1  data  for  Target  M 


-20  0  20  40 

Position  (Inches) 

(b)  S22  data  for  Target  M 


-60  -40  -20  0  20  40  60 

Position  (Inches) 

(e)  S22  data  for  Target  M 

Figure  12,  Preprocessed  GPR  data  for  Target  M  (vertical  MK23), 
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ELF  vs  Posibon  for  614  273<Jeg 


(a)  ELF  for  Target  M 


ETO  vs  Position  for  B14  27We9 


(b)  ETO  for  Target  M 


Figure  13,  Extracted  ELF  and  ETO  for  Target  M  (vertical  MK23)  as  a  function  of 
antenna  position. 
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S11  161b  shotput” 


^q-3  S21  161b  shotput’ 

,1  0 


l0.8 


l0  6 


l0  4 


02 


■60  -40 


Position  (Inches) 


(a)  Si  1  data  for  Target  G 

S22  161b  shotput’ 


-20  0  20  40 

Position  (Inches) 

(b)  S21  data  for  Target  G 


■60  -40  -20  0  20  40  60 

Position  (Inches) 


(c)  S22  data  for  Target  G 

Figure  14  Preprocessed  GPR  data  for  Target  G  (161b  shotput). 
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fELF  vs  Position  for  A1 5  273cieg 


(a)  ELF  for  Target  G 


ETO  vs  Position  for  A1 5  273deg 


(b)  ETO  for  Target  G 

Figure  15,  Extracted  ELF  and  ETO  for  Target  G  (161b  shotput)  as  a  function  of  antenna 
position. 
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(a)  Sii  data  for  Target  H 


S21  Item  H 
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(b)  S21  data  for  Target  H 


S22  Item  K 


Position  (Inches) 


(c)  S22  data  for  Target  H 


Figure  16  Preprocessed  GPR  data  for  Target  H  (horizontal  plate,  4"x  4"x  0,25"), 
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fELF  vs  PosiUon  for  A14  273<Jeg 


(a)  ELF  for  Target  H 


ETO  vs  Position  for  A14  273<Jeg 

I -  - - ^ ' - - - 
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S  50‘ 
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-SO  -AO  -30  -20  -10  6  10  20  30  40  50 

Position  (inches) 

(b)  ETO  for  Target  H 

Figure  17  Extracted  ELF  and  ETO  for  Target  H  (horizontal  plate,  4"x  4"  x  0,25")  as  a 
function  of  antenna  position. 
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Delay  Time  (ns) 


S1 1  Item  J 


-60  -40  -20  0  20  40  60 

Position  (inches) 

(a)  Sn  data  for  Target  J 


S21  Item  J 


Position  (Inches) 


(b)  S21  data  for  Target  J 


-60  -40  -20  0  20  40  60 


Position  (inches) 

(c)  S22  data  for  Target  J 

Figure  18,  Preprocessed  GPR  data  for  Target  J  (horizontal  plate,  8"x  8"  x  0,25"), 
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fELF  vs  Position  for  C1 4  273deg 


(a)  ELF  for  Target  J 


ETO  vs  Position  for  C14  273deg 


(b)  ETO  for  Target  J 


Figure  19,  Extracted  ELF  and  ETO  for  Target  J  (horizontal  plate,  8"x  8"  x  0,25")  as  a 
function  of  antenna  position. 
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Chapter  4  Feature  Extraction  and  UXO 

Classification  of  Blossom  Point  Targets 

4.1  Known  Target  Feature  Extraction 

The  linearity  (ELF),  orientation  (ETO)  and  length  (ETL)  features  for  all  blind  test  targets  are 
chosen  from  a  position  that  is  within  the  best  linearity  region,  i.e.  the  highest  ELF  values  among  all  the 
positions  and  all  the  passes  associated  with  that  particular  target.  This  occurs  unless  the  processing  has 
excluded  the  apparent  linearity  as  spurious,  as  in  the  case  of  a  compact,  near-surface  target,  which 
always  excites  the  Sn  response  more  strongly  than  822-  For  a  EIXO  that  does  not  have  a  large 
inclination  angle,  this  choice  gives  the  best  orientation  and  length  estimation  (LEN).  A  relative 
permittivity  of  20  was  used  in  converting  from  the  delay  time  to  depth,  and  from  the  resonant  frequency 
to  length.  This  relative  permittivity  value  was  found  to  be  approximately  the  average  value  of  the  actual 
probe  data  as  shown  in  Appendix  B. 

The  estimated  depth  (DEP),  target  length  (ETL)  and  target  orientation  (ETO)  for  the  known  targets 
were  found  to  agree  well  with  their  true  values  as  shown  in  Table  3.  It  was  found  that,  for  some  reason, 
target  D  and  E  seemed  to  be  swapped  based  on  the  depth  and  length  features.  It  should  be  noted  that  for 
a  non-linear  target  such  as  a  sphere,  the  estimated  length  is  not  directly  related  to  its  dimension  and 
therefore  has  a  larger  discrepancy.  Furthermore,  the  depth  estimation  for  a  UXO  that  has  a  large 
inclination  angle  is  not  entirely  straightforward,  as  for  example  the  vertical  MK23  (Known  Target  M) 
shown  in  Figure  20  (same  as  Figure  12(a)).  Figure  21  also  plots  the  corresponding  bipolar  responses  for 
comparison.  Since  the  stronger  responses  can  only  be  obtained  from  oblique  angles,  the  middle  portion 
of  the  hyperbolic  arc  is  weak  and  usually  not  obvious.  Since  the  target  depth  is  estimated  from  the 
delay  time  measured  from  the  peak  of  the  hyperbolic  arc,  extrapolation  was  used  as  indicated  by  the 
small  circles.  Notice  the  two  groups  of  arcs  that  do  not  have  the  same  curvature  as  indicated  by  the  blue 
and  red  dashed  lines.  The  lower  arc  (red  line)  that  has  a  smaller  curvature  is  caused  by  the  scattering 
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from  the  bottom  nose  of  the  MK23.  The  upper  are  (blue  line)  that  has  a  larger  curvature  is  caused  by 
the  scattering  from  the  top  tail  of  the  MK23.  One  cannot  always  see  both  scattering  mechanisms 
clearly.  This  depends  on  the  nose  and  tail  geometry  as  well  as  the  soil  absorption  and  length  of  the 
target.  The  ground  surface  is  located  at  about  the  -6  ns  position.  The  scattering  caused  by  disturbed 
soil  is  also  visible  near  the  surface.  Notice  that  the  above  discussion  has  utilized  the  early-time 
scattering  characteristics,  which  was  not  done  in  processing  for  our  previous  Demo.  Figure  13(a) 
shows  that  the  largest  ELF  values  occur  at  offset  positions  for  this  vertical  UXO.  For  other  nonlinear 
targets  or  UXO-like  targets  with  near  vertical  inclination  angles,  the  ETO  will  be  the  same  as  the  scan 
direction  for  the  reason  discussed  in  Section  3.2  and  does  not  have  much  meaning. 


Position  (inches) 


Figure  20,  Energy  plot  of  the  measured  Sll  data  for  the  vertical  MK23  (Known  Target  M), 
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S,, ,  Target  M 


Antenna  Position  (inches) 

Figure  21  Bipolar  plot  of  the  measured  Sll  data  for  the  vertical  MK23  (Known  Target  M), 


Table  3  Features  Extracted  for  Known  Targets 


TAR# 

UXO 

Confidence 

ETL 

(m) 

True  L 

ETO 

True 

Orient. 

DEP 

(m) 

True 

Dep. 

A13 

Y 

M 

0.448 

c 

12 

0 

0.246 

0.750 

A14 

N 

H 

0.196 

* 

-90 

* 

0.064 

0.050 

A15 

N 

H 

0.131 

* 

-92 

* 

0.195 

0.250 

B13 

Y 

M 

0.203 

0.150 

-84 

* 

0.137 

0.100 

B14 

Y 

H 

0.197 

* 

-94 

* 

0.202 

0.250 

B15 

Y 

H 

0.381 

0.305 

-8 

0 

0.289 

0.250 

C13 

Y 

L 

0.137 

0.305 

-92 

* 

0.195 

0.500 

C14 

N 

H 

0.239 

* 

-97 

* 

0.210 

0.250 

C15 

Y 

H 

0.152 

0.150 

0 

0 

0.253 

0.200 

D13 

Y 

H 

0.116 

* 

-87 

* 

0.340 

0.350 

D14 

Y 

L 

0.303 

0.305 

-91 

* 

0.202 

0.350 

D15 

N 

L 

0.153 

0.076 

89 

* 

0.042 

0.050 

E14 

Y 

H 

0.430 

0.610 

-9 

0 

0.631 

0.750 

E15 

Y 

L 

0.243 

0.150 

-89 

* 

0.246 

0.080 

Graphically,  we  show  Table  3  table  in  histogram  form  below  in  Figure  22.  In  the  ease  of  target 
length  only  the  eases  are  shown  for  whieh  length  determination  applies.  Cases  are  sorted  in  either 
length  or  depth  magnitude  order,  not  by  target  number  or  letter.  These  plots  illustrate  that  the  system 
for  determination  of  length  and  depth  parameters  is  basieally  sound. 


Comparison  of  known  target  lengths  to  estimated  length  (LEN) 


length 

(m) 


Comparison  of  known  target  depths  to  estimated  depth  (DEP) 


depth 

(m) 


Figure  22,  Comparisons  of  true  and  estimated  known  targets  lengths  (top)  and  depths  (bottom), 
sorted  hy  magnitude  order. 
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4.2  Offset  Tests  for  Effects  of  Positional  Uncertainty 


Additional  data  were  taken  for  the  known  target  set  with  scan  center  offset  from  the  target  location 
designated  by  NRL.  Two  different  offsets,  10  inches  and  20  inches,  were  taken  to  the  west  and  to  the 
south  with  respect  to  the  designated  target  center.  Figure  23  shows  an  example  of  10-inch  offset  for 
target  B15,  a  north-oriented  steel  strip  (3”xl2”x0.25”)  buried  at  a  depth  of  25cm  (-'10  inches).  In 
terms  of  diagonal  distance  to  the  intersection  of  perpendicular  survey  lines,  the  two  offsets  correspond 
to  about  14  and  28  inches.  The  question  we  ask  here  is  whether  the  ability  of  our  system  to  obtain  key 
processing  parameters,  discussed  above,  is  degraded  by  such  offsets.  In  our  parlance  here,  this  strip  is  a 
UXO-like  object  in  that  it  presents  an  aspect  ratio  of  at  least  3  to  1  in  all  orientations.  The  “good” 
performance  for  the  offset  survey  lines  would  mean  extracting  signal  parameters  corresponding  to  a 
UXO-like  classification.  The  extracted  classification  features  and  the  responses  for  polarization 
parallel  to  the  scan  direction  are  shown  in  Figure  24.  The  transverse  responses  are  much  weaker 
compared  to  the  parallel  ones  and  are  not  shown.  The  classification  features  clearly  indicated  a  UXO- 
like  object,  without  degradation  due  to  offset.  The  high  ELF  region  is  offset  from  the  survey  line 
midpoint,  as  expected.  The  same  target  was  measured  again  with  the  offset  increased  to  20  inches 
(equivalent  diagonal  offset  of  28  inches),  as  illustrated  in  Figure  25.  Again,  with  a  positional  shift 
corresponding  to  the  offset,  the  classification  parameters  still  clearly  show  UXO-like  features  (see 
Figure  26).  Therefore,  for  horizontal  ordnance,  it  seems  that  position  offsets  up  to  20  inches  do  not 
create  a  problem  for  classification  using  the  existing  features  and  criteria.  This  conclusion  was  also 
found  to  be  valid  for  vertical  ordnance  as  shown  in  Figure  27  and  Figure  28.  While  the  patterns  of  ELF 
and  ETO  in  these  plots  may  appear  a  bit  ragged  for  those  not  accustomed  to  the  typical  variability  of 
these  signals,  the  general  patterns  contrast  notably  to  cases  lacking  significant  linearity,  or  failing  to 
show  coincidences  of  linearity  and  distinct  orientation.  Note  that  a  given  offset  should  cause  less  of  a 
problem  for  deeper  targets,  assuming  that  increased  signal  absorption  over  depth  does  not  obscure  the 
data  content.  Overall,  discounting  attenuation  effects  over  depth,  we  expect  offset  effects  to  be  similar 
for  cases  when  the  ratio  of  offset  to  depth  is  similar.  In  the  case  of  target  B 1 5,  good  results  were 
obtained  for  offset  (closest  point  on  survey  path  to  target  location)equal  to  twice  the  target  depth. 
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Figure  23,  A  north-directed  scan  with  its  scan  center  offset  from  target  center  hy  10  inches  from 
both  east  and  north  direction. 


fELF  vs  Position  for  B15 


ETO  vs  Position  for  B1 5  ATV  =  357 


(a)  Estimated  linearity  faetor  (ELF)  (b)  Estimated  target  orientation  (ETO)  relative  to  North 
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Figure  24 


Classification  features  and  position-time  plot  for  target  B15,  obtained  form  the  10- 
inch  (14  inch  diagonal)  offset  survey  path. 
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Figure  25  An  East-West  scan  with  its  scan  center  offset  from  target  center  by  20  inches  from 
both  east  and  north  direction. 


fELF  vs  Position  for  B1 5  ETO  vs  Position  for  B 1 5  ATV  =  268 


(a)  Estimated  linearity  function  ELF  for  (b)  Estimated  target  orientation  ETO  for  target 

target  B 1 5  with  survey  path  offset  20  inches.  B 1 5  with  survey  path  offset  20  inches. 


C  -30 


-60  -40  -20  0  20  40 

Position  (inches) 


S22  B1 5  Orientation  =  268 


(c)  Sii  (transmitted  and 
received  signals  both 
parallel  to  scan  direction), 
as  a  function  of  antenna 
position. 


Figure  26  Classification  features  and  position-time  signal  plot,  for  target  B15,  obtained  from  the 
20-inch  offset  survey  path. 
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fELF  vs  Position  for  B1 4 


ETO  vs  Position  for  B14  ATV  =  181 


(a)  Estimated  linearity  function  ELF  for  vertical 
MK23  with  survey  path  offset  10  inches. 


(b)  Estimated  target  orientation  ETO  for  vertical 
MK23  with  survey  path  offset  10  inches 


S11  B14  Orientation  =  181 


Position  (inches) 


(c)  Sii  (transmitted  and  received  signals  both  parallel  to  scan  direction),  as  a  function  of  antenna 
position. 


Figure  27  Classification  features  for  a  vertical  MK23  bomb  obtained  form  tbe  10-incb  offset 
configuration. 
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fELF  vs  Position  for  BM 


(a)  Estimated  linearity  function  ELF  for  vertical 
MK23  with  survey  path  offset  20  inches. 


ETO  vs  Position  for  B 1 4  ATV  =  1 0 


(b)  Estimated  target  orientation  ETO  for  vertical 
MK23  with  survey  path  offset  10  inches 


S11  B14  Orientation  =  10 


-20  0 
Position  (inches) 


(c)  Sii  (transmitted  and  received  signals  both  parallel  to  scan  direction),  as  a  function  of  antenna 
position,  for  20  inch  offset. 


Figure  28  Classification  features  and  position-time  signal  plot  for  a  vertical  MK23  bomb 
obtained  from  a  20-incb  offset  (28  inch  diagonal)  survey  path. 
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4.3  UXO  Classification  Algorithm 


The  previous  UXO  elassifieation  algorithm  utilized  resonanee  and  linearity  features  based  on  single 
position  GPR  data.  Under  sueh  a  eriteria,  the  offset  of  shallow  nonlinear  objects  caused  false  alarms. 
UXO’s  with  large  inclination  angles  were  easily  missed.  This  situation  has  been  improved  significantly 
after  the  introduction  of  spatial  and  scattering  features  obtained  from  multiple  position  data.  Figure  29 
shows  the  flow  chart  of  the  new  UXO  classification  criteria. 

The  whole  UXO  classification  procedure  starts  with  an  inspection  of  the  spatial  distributions  of  the 
extracted  ELF,  that  is,  of  the  ELF  plotted  as  a  function  of  antenna  positions.  If  the  ELF  is  low  over 
most  of  the  scan  region  (10  feet),  it  means  the  object  is  very  likely  a  non-UXO  object.  If  there  are  some 
regions  of  high  ELF  values,  the  next  thing  checked  is  whether  high  ELF  values  are  observed  in  the 
middle  region  such  as  the  one  shown  in  Figure  1 1(a).  If  the  answer  is  yes,  then  the  object  is  very  likely 
to  be  a  UXO-like  object  with  a  small  inclination  angle.  The  corresponding  ETO  indicates  the  UXO 
orientation.  If  the  ELE  values  in  the  middle  region  are  low  but  high  at  offset  regions  (whether  one  side 
or  two  sides),  the  object  could  be  a  UXO  with  a  larger  inclination  as  shown  in  Eigure  13(a)  or  a  shallow 
non-linear,  non-UXO  object  as  shown  in  Eigure  15(a),  Eigure  17  and  Eigure  19.  The  former  could  have 
a  higher  ELE  due  to  oblique  incidence.  The  latter  produces  a  higher  ELE  because  of  coupling  with  the 
antenna  arms.  In  either  situation,  ETO  will  indicate  an  orientation  aligned  with  the  scan  direction  in  all 
passes  as  seen  in  Eigure  13(b)  and  Eigure  15(b)  where  the  ETO’s  in  both  cases  indicate  an  angle  of- 
90°  that  agrees  with  the  scan  direction,  273°. 

A  plate-like  object,  such  as  shrapnel,  may  show  good  resonance  depending  on  its  geometry  and 
edges.  However,  the  surface  reflection  from  a  horizontal  plate-like  object  produces  strong  responses  in 
both  Sii  and  S22  channels,  as  demonstrated  in  Eigure  18  above.  On  the  other  hand,  a  vertical  UXO 
would  have  much  weaker  S22  responses  as  shown  in  Eigure  12  because  the  field  polarization  is 
transverse  to  the  UXO  axis.  This  quite  different  S22  channel  behavior  allows  one  to  discriminate  a 
plate-like  object  from  a  vertical  one  without  too  much  difficulty.  It  was  also  observed  that  the  scattering 
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pattern  in  time-position  plot  under  the  initial  ares  also  shows  a  more  complieated  structure  compared  to 
a  vertical  UXO.  This  could  be  related  to  the  different  creeping  waves  propagating  behind  the  plate. 


A  compact  non-linear  object  such  as  a  sphere  is  not  a  good  resonator  and  can  be  discriminated  by 
the  large  damping  factor  extracted  from  the  natural  resonance.  This  type  of  object  also  does  not  have 
the  two  arcs  with  different  curvature  as  observed  in  Figure  20. 

The  classification  criteria  discussed  above  were  found  to  be  very  effective  in  discriminating  UXO- 
like  targets  for  the  lanes  with  known  targets.  Each  of  these  criteria  may  be  developed  into  automatic 
classification  procedures  using  pattern  recognition,  image  correlation  or  neural  network  training 
techniques.  However,  at  present.  Figure  29  is  implemented  by  training  an  operator  using  a  training  set. 
The  operator  is  then  asked  to  make  a  classification  decision  by  following  the  flow  chart. 
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Object  Large  Inclination  Angle 


Figure  29,  Improved  UXO  classification  fiow  chart. 


39 


Chapter  5  Blind  Targets  Classification  Results 
Compared  with  Ground  Truth 


This  section  starts  with  the  designation  of  UXO-like  (or  ordnance)  and  clutter  (or  non-ordnance) 
items.  The  overall  classification  performance  will  then  be  provided  for  both  new  and  old  lanes  at  the 
Blossom  Point  site.  The  whole  target  list  is  then  broken  down  into  three  parts  for  further  investigation: 
correctly  classified  UXO-like  items,  missed  UXO-like  items  and  clutter  items  classified  as  UXO-like 
items. 


5.1  Designation  of  UXO-Like  Items 

Before  we  evaluate  the  classification  performance,  all  the  target  items  will  be  separated  into  UXO- 
like  and  clutter  two  categories.  It  should  be  noted  that  some  non-ordnance  items  may  have  general 
elongated  bodies  like  a  typical  UXO  item  regardless  their  different  detail  shapes,  smoothness  or 
curvature.  Conductive  cylinders  or  shrapnel  strips  are  examples.  This  type  of  non-ordnance  would 
show  similar  radar  features  under  the  practical  frequency  limitations  and  current  classification  criteria. 
Therefore,  when  an  item  has  a  length-to-diameter  (L/D)  ratio  of  three  or  greater,  it  is  designated  as  an 
UXO-like  item. 

Table  4  lists  the  UXO-like/Clutter  designation  of  the  old  Blossom  Point  items  buried  in  lanes  A-E, 
and  Pictures  of  these  items  can  also  be  found  in  Figure  30.  Some  "new"  items  were  added  to  the  site 
specifically  for  this  demo,  including  standard  ordnance  and  various  odd  pieces  of  scrap  from  exploded 
ordnance.  These  1 1  scrap  items  are  numbered  here  Scl  through  Sell.  Some  of  these  bear  cataloguing 
numbers  from  previous  work,  indicated  on  blue  tags  attached  to  each  piece.  One  target  (Sc4)  consists 
of  a  cluster  of  pieces  with  tag  numbers  12,  17,  and  25.  Vertical  and  horizontal  views  of  these  pieces  are 
shown  in  Figure  31.  Table  5  shows  the  UXO-like/clutter  designation  of  the  new  Blossom  Point  items 
buried  in  the  F-'I  lanes. 
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Table  4,  UXO-Like/Clutter  Designation  of  Old  Blossom  Point  Items  (Lanes  A~E), 


Item  ID 

Description 

Classification 

Reason 

A 

1 .5"x3"  solid  steel  cylinder 

Clutter 

L/D<3 

B 

1 .5"x  6"  solid  steel  cylinder 

UXO-Like 

L/D>  3 

C 

1.5"x12"  solid  steel  cylinder 

UXO-Like 

L/D>3 

D 

3"x6"  solid  steel  cylinder 

Clutter 

L/D<3 

E 

3"  X  12"  solid  steel  cylinder 

UXO-Like 

L/D>  3 

F 

3"  X  24"  solid  steel  cylinder 

UXO-Like 

L/D>3 

G 

16  lb  steel  shotput 

Clutter 

Symmetric 

H 

4"x4"x0.25"  steel  plate 

Clutter 

Symmetric 

J 

8"x8"x0.25"  steel  plate 

Clutter 

Symmetric 

K 

1 .5"x6"x0.25"  steel  plate 

UXO-Like 

L/D>3 

L 

3"x12"x0.25"  steel  plate 

UXO-Like 

L/D>  3 

M 

Mk  23  practice  bomb 

UXO-Like 

UXO 

N 

8 1  mm  Mortar 

UXO-Like 

UXO 

P 

BDU-33 

UXO-Like 

UXO 

Qi 

4"x4"x0.25"  A1  plate 

Clutter 

Symmetric 

Q2 

2"x8"x0.25"  A1  plate 

UXO-Like 

L/D>3 

Q3 

Flattened  A1  soda  can 

Clutter 

L/D<3 

Q4 

Box  fin  (M  38) 

Clutter 

Not  Linear 

Q5 

Box  fin  (M  38) 

Clutter 

Not  Linear 

Q6 

Slip  joint  pliers 

UXO-Like 

L/D>3 

Q7 

A1  fins  (2.75"  rocket) 

UXO-Like 

L/D>3 

Q8 

Shovel  blade 

UXO-Like 

L/D  >  3  (handle) 

Q9 

Banding  material 

Clutter 

L/D<3 

QIO 

Box  fin  assembly  (M  38) 

Clutter 

Not  Linear 

Qll 

Barbed  wire 

Clutter 

Not  Linear 

Q12 

Florseshoe 

Clutter 

Not  Linear 

Q13 

8”  round  steel  plate 

Clutter 

Symmetric 

Q14 

Twisted  steel  piece 

Clutter 

L/D  <3 

Q15 

Fence  post  cap 

Clutter 

L/D  <3 
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(a) 


(b) 

Figure  30,  Pictures  of  the  old  Blossom  Point  items  in  lanes  A~E,  (a)  calibration  targets, 
clutter  items 


(b) 
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Table  5,  UXO-Like/Clutter  Designation  of  New  Blossom  Point  Items  (Lanes  F~I) 


Item  ID 

Description 

Designation 

Reason 

X1 

60mm 

UXO-Like 

UXO 

X2 

60mm 

UXO-Like 

UXO 

X3 

60mm 

UXO-Like 

UXO 

X4 

60mm 

UXO-Like 

UXO 

X5 

60mm 

UXO-Like 

UXO 

X6 

2.75" 

UXO-Like 

UXO 

X7 

2.75" 

UXO-Like 

UXO 

X8 

2.75" 

UXO-Like 

UXO 

X9 

2.75" 

UXO-Like 

UXO 

X10 

81mm 

UXO-Like 

UXO 

X11 

81mm 

UXO-Like 

UXO 

X12 

81mm 

UXO-Like 

UXO 

X13 

BDU-33 

UXO-Like 

UXO 

X14 

120mm 

UXO-Like 

UXO 

NX1 

Steel  Bar  (L24"/D3") 

UXO-Like 

L/D>3 

NX2 

Steel  Bar  (UXO  Like  L=same  as  120mm, 
D1" 

UXO-Like 

L/D>3 

NX3 

8"d  steel  plate 

Clutter 

Symmetric 

NX4 

8"d  steel  plate 

Clutter 

Symmetric 

NX5 

8"d  steel  plate 

Clutter 

Symmetric 

Sc1 

Whole  crumpled  box  fin 

Clutter 

Not  Linear 

Sc2 

box  fin  fragment 

Clutter 

Not  Linear 

Sc3 

Projectile  nose 

Clutter 

Not  Linear 

Sc4 

Cluster  of  clutter  items  25,17,  &12 

Clutter 

Not  Linear 

Sc5 

clutter  Item  #22 

Clutter 

Not  Linear 

Sc6 

Clutter  Item  #3 

UXO-Like 

L/D>3 

Sc7 

clutter  Item  #24 

Clutter 

Not  Linear 

Sc8 

clutter  Item  #18 

Clutter 

Not  Linear 

Sc9 

clutter  Item  #23 

Clutter 

Not  Linear 

Sc10 

clutter  Item  #1 1 

UXO-Like 

L/D>3 

Sc11 

clutter  Item  #16 

UXO-Like 

L/D>3 
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Sc6,  Sc  10,  and  Sell  are  all  clutter  items  that  are  designated  as  UXO-like  beeause  of  their  L/D  ratio 
greater  than  three.  Pictures  of  each  item  numbered  Scl  thru  Sell  can  be  seen  below. 
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Figure  31  Pictures  of  the  new  Blossom  Point  items  in  lanes  F~I, 
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5.2  Overall  Classification  Performance 


The  following  three  parameters  are  used  to  evaluate  the  overall  elassifieation  performanee. 

The  True-UXO-Rate  (TUR)  is  defined  as 

TUR  =  [Number  of  UXO-like  items  classified  as  UXO-like]  /  [Number  of  UXO-like  items] 

The  False-UXO-Rate  (FUR)  is  defined  as 

FUR  =  [Number  of  clutter  items  classified  as  UXO-.like]  /  [Number  of  clutter  items] 

The  True-Clutter-Rate  (TCR)  is  defined  as 

TCR  =  [Number  of  clutter  items  classified  as  clutter]  /  [Number  of  clutter  items] 

Notice  that  FUR+TCR=100%  and  TCR  is  directly  related  to  the  cost  reduction  for  the  clearance  task. 
The  risk  factor,  or  the  UXO-like  that  are  not  classified  correctly  can  be  obtained  from  (1-TUR). 

Figure  32  plots  the  TUR,  FCR  and  TCR  for  three  different  thresholds  based  on  the  confidence  levels 
(see  Table  6)  for  all  the  targets  measured  at  the  Blossom  Point  site.  These  targets  include  the  old  items 
buried  in  lanes  A  thru  E  and  the  new  items  selected  by  CRREL/  OSU  team  buried  in  lanes  F  thru  I.  The 
leftmost  group  in  the  figure  neglects  of  the  confidence  levels  and  is  solely  determined  by  whether  the 
item  is  classified  as  UXO-like  or  not.  The  middle  group  includes  only  those  UXO-like  items  that  have 
moderate  or  high  confidence  levels  as  UXO-like  items.  All  other  items  are  classified  as  clutter  items. 
The  rightmost  group  includes  only  those  UXO-like  items  that  have  high  confidence  levels  as  UXO-like 
items.  Although  more  than  75%  of  UXO-like  items  were  correctly  classified  as  indicated  by  the  (TUR), 
further  improvement  is  needed.  The  number  of  clutter  items  that  were  classified  as  UXO-like  is 
approximately  48%  as  indicated  by  FUR.  Both  TUR  and  FUR  can  be  further  improved  as  will  be 
address  in  next  section. 


47 


It  was  apparent  during  the  data  proeessing  that  the  lanes  eontaining  the  new  targets  suffered  from 
mueh  higher  surfaee  clutter  caused  by  recent  soil  disturbance.  It  was  also  apparent  during  the 
demonstration  that  there  were  many  parallel  linear  trenches  running  in  the  north-south  and  east-west 
directions.  The  trenches  oriented  in  the  north-south  direction  have  a  width  of  approximately  14  inches. 
The  east-west  trenches  are  narrower.  These  trenches  could  be  caused  by  the  extraction  of  previous 
drainage  or  irrigation  pipes  installed  while  the  site  was  still  a  farmland.  The  new  items  also  contain 
more  target  orientation  variation.  With  these  factors  in  mind,  it  is  revealing  to  examine  classification 
performance  for  old  items  and  new  items  separately.  Figure  33  plots  the  TUR,  FUR  and  TCR  for  old 
Blossom  Point  items  buried  in  lanes  A  thru  E.  All  of  these  items  are  either  horizontal  or  vertical.  The 
classification  results  without  using  any  threshold  shows  amazingly  that  only  one  UXO-like  item  is 
missed!  The  FUR  is  lower  compared  to  the  previous  case  but  is  still  45%.  The  classification  results  for 
the  new  target  items  alone  are  shown  in  Figure  34,  indicating  a  much  poorer  performance.  It  should  be 
noted  that  there  are  39  UXO-like  items  and  22  clutter  items  in  the  old  lanes.  There  are  19  UXO-like 
items  and  1 1  clutter  items  in  the  new  lanes.  The  following  sections  will  investigate  further  the  issues  of 
high  FUR  and  missed  ordnance.  These  results  include  the  known  targets  located  at  Rows  n-lS  in 
Lanes  A~E.  During  the  classification,  the  operator  did  not  pay  any  attention  to  whether  a  target  was 
known  or  unknown  and  performed  classification  solely  based  on  classification  features.  While 
inclusion  of  the  calibration  targets  in  the  scoring  may  be  regarded  as  suspect,  note  in  Eigure  33  and 
Eigure  36  that  the  histograms  for  the  old  lanes  are  similar  whether  or  not  the  calibration  targets  are 
included.  Note  that  the  known  targets  are  only  located  in  the  old  lanes.  Therefore  the  statistics  for  new 
lanes  are  not  affected  by  inclusion  or  deletion  of  data  from  the  cal  lanes. 
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Entire  MTR  Site 


None  Low  Moderate 

Threshold 


Figure  32,  The  classification  performance  for  all  targets  buried  at  the  Blossom  Point  site 
considered  together  (old  ones  lanes  A~E  and  new  ones  lanes  F~I),  including 
calibration  targets  treated  as  unknowns. 


Original  MTR  Items 


None  Low  Moderate 


Threshold 


Figure  33  Classification  performance  for  the  old  targets  alone  (A~E  lanes)  buried  at  the  Blossom 
Point  site,  including  calibration  targets  treated  as  unknowns. 
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New  GPR  Items 


100 


None  Low  Moderate 

Threshold 

Figure  34,  The  classification  performance  for  the  new  targets  (F~I  lanes)  hurled  at  the  Blossom 
Point  site, 
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Figure  35,  The  classification  performance  for  all  lanes  excluding  the  calibration  targets. 
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Original  MTR  Items 

100^ - 


None  Low  Moderate 

Threshold 


Figure  36,  The  classification  performance  for  the  old  lanes  excluding  the  calibration  targets. 


5.3  Correctly  Classified  UXO-like  Items 

Table  6  lists  the  correetly  classified  UXOs,  i.e.  TUR  items  in  comparison  with  the  ground  truth 
information  provided  by  NRL.  The  last  column  indicates  the  signal-to-clutter  ratio  SCR  obtained  from 
the  measured  data.  The  signal  level  of  the  SCR  was  selected  from  the  maximum  eigenvalues  associated 
with  the  selected  target  among  different  channels  and  positions.  The  clutter  level  was  determined  from 
the  eigenvalues  obtained  from  the  empty-site  data  using  exact  the  same  processing  procedures  and 
parameters  used  in  obtaining  the  signal  level. 

First,  let’s  compare  the  estimated  depths  with  their  true  values  as  shown  in  Figure  37.  A  general 
agreement  for  most  targets  is  observed.  Items  A13,  B7,  Cl 3,  D4,  Dll  and  18  are  the  ones  that  have 
greater  depth  discrepancy.  These  are  also  highlighted  in  Table  6.  Notice  the  large  discrepancies 
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observed  for  items  18  and  D4.  Both  items  are  vertical  UXO-like  items  buried  relatively  deep  (one 
meter).  The  combination  of  weak  scattering  from  this  orientation  as  well  as  the  soil  absorption  caused 
the  responses  to  be  weak  compared  to  the  background  clutter  level.  From  the  estimated  depths  (DEP),  it 
is  obvious  that  the  processing  actually  picked  up  surface  anomalies  in  preference  to  the  target 
responses.  At  this  particular  target  location,  there  was  also  a  linear  trench  running  approximately  south- 
north,  that  was  confirmed  by  the  GPR  data  and  a  log  taken  at  the  site.  This  type  of  trench  will  give  a 
stronger  S22  response  and  result  in  a  linear  feature  with  an  orientation  perpendicular  to  the  trench 
direction.  This  is  confirmed  by  the  estimated  orientations  (ETO)  obtained  from  these  two  items.  It  is 
also  interesting  to  notice  that  all  of  these  targets  have  either  low  or  moderate  confidence  levels  and 
relatively  low  signal-to-clutter  ratio  (SCR).  Thus,  while  erroneous  classification  results  were  obtained, 
the  system  indicated  that  the  results  should  be  regarded  with  suspicion.  Eurther,  while  the  processing 
missed  targets  that  were  present,  it  was  correct  in  the  sense  that  it  classified  signals  from  non-targets  as 
non-UXO.  This  motivates  us  to  believe  that  the  results  can  be  improved  by  integration  of  non-GPR 
data  on  depth  estimations,  so  that  selection  of  appropriate  time  regions  is  more  likely  to  occur  in  the 
processing.  In  the  least,  it  will  be  worth  it  to  segregate  cases  such  that  deeper  and  shallower  signals  can 
be  treated  separately. 

The  estimated  lengths  are  compared  with  true  lengths  in  Eigure  38.  Recall  that  length  estimation  was 
obtained  from  the  natural  resonance  frequency.  The  error  sources  of  length  estimation  include  incorrect 
resonant  frequency  estimation  and  soil  property  variation.  The  error  in  natural  resonant  frequency 
estimation  could  result  from  low  signal  to  clutter  ratio  and  selection  of  inappropriate  time  regions  as 
“late-time.”  The  length  estimations  for  most  selected  UXO-like  items  are  close  to  their  true  values. 

The  items  that  have  greater  length  discrepancy  are  B7,  B9,  Cl 3,  D4,  D1 1,  D13,  E2,  E14  and  113  as  also 
highlighted  in  Table  6.  It  is  not  surprising  that  B7,  CI3,  D4  and  DI I  have  larger  length  estimation 
errors  as  their  depths  were  estimated  incorrectly  and  they  all  have  low  confidence  levels.  Item  DI3  and 
EI4  have  similar  lengths  (-O.bm)  and  seem  to  have  similar  underestimation  of  the  lengths.  The  lengths 
of  items  B7,  Cl 3,  D4  and  D13  are  all  greatly  underestimated.  This  could  be  explained  by  the  fact  that 
they  all  had  vertical  orientations  and  the  higher  resonant  modes  were  excited  more  strongly  than  the 
fundamental  mode  (half-wavelength  mode  used  for  length  estimation)  due  to  the  oblique  incidence. 

This  resulted  in  overestimation  of  the  resonant  frequency,  which  in  turn  caused  the  target  length  to  be 
underestimated.  Erom  examination  of  the  scattering  patterns  over  space,  it  appears  that  such  an 
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underestimation  is  more  likely  for  longer  objeets,  as  the  ineident  angle  gets  eloser  to  the  ordnanee  axis 
direetion.  Therefore,  it  may  not  be  praetieal  to  estimate  the  length  of  vertieal  ordnanee  aeeurately  unless 
it  is  shallow  and  the  soil  is  not  too  absorptive. 

Figure  39  eompares  the  estimated  target  orientations  with  the  true  orientations,  for  orientable  targets. 
Sinee  the  orientation  of  a  vertieal  item  does  not  have  mueh  meaning,  all  vertieal  items  have  been 
exeluded  from  this  figure.  An  exeellent  orientation  agreement  is  observed.  Only  item  C8,  D3. 15  and 
113  show  error  greater  than  15°. 


Table  6  List  of  correctly  classified  UXO-like  items. 


Confidence 

Level 

Est. 

True 

Est. 

True 

Est. 

True 

True 

TAR# 

Length 

(m) 

Length 

(m) 

Depth 

(m) 

Depth 

(m) 

Angle 

(+/-180") 

Azimut 

hO 

Dipping 

(") 

SCR 

A1 

H 

0.536 

0.51 

0.297 

0.35 

-174 

0 

0 

48.46 

A2 

H 

0.371 

0.25' 

0.231 

0.23 

-177 

0 

0 

64.40 

A3 

L 

0.165 

0.305 

0.326 

0.5 

93 

0 

90 

23.82 

A4 

H 

0.540 

0.51 

0.398 

0.5 

-179 

0 

0 

54.95 

A5 

H 

0.276 

0.305 

0.348 

0.35 

180 

0 

0 

39.52 

A7 

H 

0.258 

0.305 

0.464 

0.5 

-176 

0 

0 

30.52 

A12 

M 

0.164 

0.152 

0.297 

0.1 

84 

0 

90 

48.99 

A13 

M 

0.448 

0.6 

0.246 

0.75 

12 

0 

90 

21.80 

B2 

M 

0.277 

0.254 

0.529 

0.15 

-173 

0 

0 

43.57 

B4 

H 

0.346 

0.305 

0.326 

0.35 

18 

0 

90 

34.28 

B7 

L 

0.255 

0.61 

0.246 

0.75 

-163 

0 

90 

20.96 

B9 

H 

0.750 

0.305 

0.260 

0.35 

176 

0 

0 

65.17 

BIO 

H 

0.293 

0.254 

0.268 

0.25 

-175 

0 

0 

35.12 

B12 

H 

0.280 

0.305 

0.180 

0.15 

178 

0 

0 

48.93 

B14 

H 

0.197 

0.254 

0.202 

0.25 

-94 

0 

90 

29.01 

B15 

H 

0.381 

0.305 

0.289 

0.25 

-8 

0 

0 

43.52 

Cl 

H 

0.190 

0.305 

0.406 

0.5 

11 

0 

0 

36.75 

C4 

H 

0.220 

0.152 

0.239 

0.2 

11 

0 

90 

21.00 

C6 

H 

0.183 

0.152 

0.151 

0.1 

7 

0 

0 

38.82 

C8 

L 

0.193 

0.152 

0.101 

0.05 

58 

90 

0 

47.84 

C13 

L 

0.137 

0.305 

0.195 

0.5 

-92 

0 

90 

32.10 

CIS 

H 

0.152 

0.152 

0.253 

0.2 

0 

0 

0 

43.33 

D2 

H 

0.302 

0.305 

0.246 

0.15 

96 

0 

90 

53.02 

D3 

H 

0.325 

0.305 

0.231 

0.25 

-132 

90 

0 

66.09 

D4 

L 

0.174 

0.61 

0.021 

1 

95 

0 

90 

17.03 

D5 

H 

0.263 

0.254 

0.202 

0.15 

9 

0 

90 

36.90 

D7 

H 

0.350 

0.305 

0.326 

0.25 

43 

0 

90 

79.78 

'  Handle  part  only. 
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Dll 

L 

0,303 

0,61 

0,667 

1 

-175 

0 

0 

31.45 

D12 

H 

0.292 

0.25 

0.115 

0.1 

-177 

0 

0 

73.96 

D13 

H 

0,116 

0,51 

0.340 

0.35 

-87 

0 

90 

32.61 

D14 

L 

0.303 

0.305 

0.202 

0.35 

-91 

0 

90 

22.05 

E2 

H 

0,363 

0,6 

0.529 

0.75 

-175 

0 

0 

45.50 

E3 

H 

0.185 

0.203 

0.130 

0.05 

-171 

0 

0 

68.51 

E9 

H 

0.185 

0.15 

0.159 

0.05 

-172 

0 

0 

53.05 

ElO 

H 

0.389 

0.305 

0.188 

0.15 

-127 

90 

0 

53.18 

E12 

H 

0.188 

0.152 

0.108 

0.05 

179 

0 

0 

44.16 

E14 

H 

0,153 

0,61 

0.631 

0.75 

-9 

0 

0 

50.55 

E15 

L 

0.243 

0.152 

0.246 

0.08 

-89 

0 

90 

17.88 

*14 

H 

0.295 

0.23 

0.420 

0.25 

175 

0 

0 

69.40 

*15 

M 

0.206 

0.23 

0.515 

0.5 

29 

45 

30 

19.30 

*18 

M 

0.210 

0.23 

0,057 

1 

81 

0 

90 

22.95 

*19 

M 

0.160 

0.42 

0.311 

0.25 

5 

0 

30 

23.76 

*110 

M 

0.307 

0.42 

0.108 

0.5 

13 

0 

45 

15.42 

*113 

H 

0,364 

0,51 

0.355 

0.25 

23 

0 

30 

50.95 

*F  to  I  lanes  are  new  sites  containing  new  GPR  targets. 
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Figure  37,  Comparison  of  estimated  depth  and  true  depth  of  correctly  classified  UXO-like  items. 
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Figure  38,  Comparison  of  estimated  length  and  true  length  of  correctly  classified  UXO-like 
items. 


Figure  39,  Comparison  of  estimated  orientation  and  true  orientation  of  correctly  classified 
UXO-like  items 
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5.4  Missed  UXO-like  Items 


Table  7  lists  the  missed  UXO-items  to  be  analyzed  in  more  detail.  It  is  interesting  to  note  that  only 
one  vertical  UXO  (B3)  in  old  lanes  was  missed  out  of  39  UXO-like  items  and  22  clutter  items.  Also 
notice  that  the  estimated  depths  for  most  of  the  missed  UXO-like  items  are  very  close  to  surface.  This 
is  probably  due  to  the  weakness  of  the  true  target  response  so  that  the  surface  anomalies  were 
mistakenly  selected  as  the  target  responses.  This  certainly  is  the  case  for  F6,  F7,  G7, 16, 17  and  112 
(highlighted  in  Table  7)  since  they  were  all  small  items  and  were  buried  relatively  deep.  The  responses 
of  items  H4  and  H5  were  virtually  indistinguishable  from  the  background  clutter. 

B3,  II 1  and  II4  were  found  to  show  UXO-like  features  and  should  be  reclassified  as  UXO-like 
items  after  choosing  the  correct  time  range  corresponding  to  the  true  depth.  The  data  collected  from  the 
new  lanes  at  the  Blossom  Point  site  showed  significantly  stronger  ground  surface  scattering  compared 
to  those  collected  at  the  old  lanes.  For  several  weak  targets,  the  predicted  depth  clearly  indicated  that 
we  have  selected  the  surface  clutter  instead  of  the  targets  responses.  In  a  sense,  this  is  a  “good” 
because  we  have  classified  this  surface  scattering  as  non-UXO! 


Table  7  List  of  missed  UXO-like  items. 


TAR 

# 

Confidence 

Level 

Est. 

Depth 

(m) 

True 

Depth 

(m) 

Est. 

Orient. 

(+/-180") 

True 

Azimut 

hC) 

True 

Dipping 

(") 

True 

Length 

(m) 

Ground  Truth  ID 

SCR 

B3 

H 

0.035 

0.5 

92 

0 

90 

0.66 

N0305004 

-8.69 

*111 

H 

0.071 

0.5 

-85 

90 

60 

0.42 

2.75" 

12.24 

*114 

H 

0.028 

0.5 

4 

45 

45 

0.51 

81mm 

12.86 

*F6 

M 

0.057 

0.5 

175 

45 

60 

0.3 

clutter  Item  #  1 1  (Sc  1 0) 

8.14 

*F7 

H 

0.079 

0.5 

5 

45 

45 

0.4 

clutter  Item  #16(Scl  1) 

-5.58 

*G7 

H 

0.064 

0.5 

85 

0 

90 

0.33 

clutter  Item  #3  (Sc6) 

-7.14 

*16 

H 

0.064 

0.75 

153 

90 

30 

0.24 

60mm 

25.66 

*17 

H 

0.064 

0.75 

-1 

90 

45 

0.24 

60mm 

-15.45 

*112 

M 

0.057 

0.75 

171 

45 

75 

0.42 

2.75" 

5.32 

*H4 

H 

0.057 

1 

127 

45 

30 

0.69 

BDU-33 

0.65 

*H5 

M 

0.159 

1 

84 

45 

60 

0.75 

120mm 

10.33 

*H6 

H 

0.086 

0.75 

42 

0 

30 

0.61 

Steel  Bar  (L24"/D3") 

9.35 

*H7 

H 

0.086 

0.75 

82 

90 

45 

0.75 

Bar  L=same  as  120mm 

23.21 

*115 

H 

0.064 

0.75 

-6 

90 

60 

0.51 

81mm 

11.57 

*  New  lanes. 
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5.5  Clutter  Classified  as  UXO-like  Items 


A  list  of  clutter  items  that  were  classified  as  UXO-like  objects  is  provided  in  Table  8.  Further 
investigations  of  these  cases  is  summarized  below. 


Table  8  List  of  clutter  classified  as  UXO-like  items 


TAR# 

Confidence 

Level 

Est. 

Depth 

(m) 

True 

Depth 

(m) 

Est. 

Orient. 

(+/-180“) 

True 

Azimuth 

(“) 

True 

Dipping 

(“) 

Ground  Truth  ID 

E8 

H 

0.035 

0.08 

-173 

0 

0 

Q1500800(Fence  Post  Cap) 

C2 

H 

0.057 

0.10 

-170 

0 

0 

Q 1 1 0 1 000  (Barbed  Wire) 

Cll 

M 

0.457 

0.10 

7 

0 

0 

Q0901000(Banding  Material) 

A9 

H 

0.166 

0.25 

-172 

0 

90 

8”x8”xl/4”  Plate 

*H11 

M 

0.042 

0.25 

-151 

Arbitrary 

Arbitrary 

crumpled  box  fin  (Scl) 

C9 

M 

0.420 

0.05 

9 

0 

0 

Horse  Shoe 

D1 

H 

0.558 

0.10 

119 

0 

0 

O0501000(Box  Fin  Piece,M38) 

D8 

H 

0.558 

0.15 

103 

0 

0 

Q1001500(Box  Fin  M38) 

*H8 

H 

0.144 

0.30 

7 

Arbitrary 

0 

8''d  steel  plate 

*H12 

L 

0.042 

0.50 

-1 

0 

0 

box  fin  fragment  (Sc2) 

*H9 

L 

0.050 

0.50 

89 

Arbitrary 

45 

8"d  steel  plate 

*F5 

L 

0.086 

0.75 

6 

90 

0 

Bomb  Shrapnel  (Sc9) 

*G9 

L 

0.057 

0.50 

177 

90 

0 

Bomb  Shrapnel  (Sc8) 

B13 

M 

0.137 

0.10 

-84 

0 

90 

3”x6”  Solid  Cylinder 

A6 

M 

0.071 

0.03 

-173 

0 

0 

Side  Crushed  Coke  Can 

CIO 

M 

0.079 

0.05 

11 

0 

90 

4”x4”xl/4”  Plate 

New  lanes. 


E8,  C2  and  Cll  -  UXO-like  geometries.  Figure  42  shows  the  classification  features  for  Item  E8.  In  this 
and  the  following  figure,  “transverse”  indicates  perpendicular  to  the  survey  pass  direction,  and 
“parallel”  indicates  aligned  with  it.  All  the  features  indicate  an  UXO-like  item  based  on  current 
classification  rules.  Although  C2  (barbed  wire)  and  Cll  (banding  material)  do  not  appear  to  have  a  L/D 
ratio  of  three  or  greater  from  the  picture,  each  individual  thread  does  have  a  thin  conducting  body. 

Since  they  are  bundled  such  that  all  threads  have  similar  orientations,  the  radar  responses  show  very 
good  linear  and  resonant  signatures. 

A9  &  HI  1  -  New  Classification  Features 
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Figure  43  shows  the  classification  features  for  Item  A9,  an  8”x8”  vertical  plate,  when  the  scan  direction 
is  transverse  to  the  top  edge  orientation.  High  ELF  values  were  observed  near  the  center  location. 
However,  unlike  a  true  UXO-like  item,  the  S 1 1  responses  are  very  weak  near  the  center  location.  This 
feature  apparently  was  not  included  in  the  existing  classification  rules  shown  in  Figure  29.  It  should  be 
noted  that  there  was  no  vertical  plate  in  the  calibration  lanes  at  the  Blossom  Point  site.  Near  the  center 
location  the  polarization  of  SI  1  is  orthogonal  to  the  edge  and  thus  produces  much  weaker  response  than 
S22  channel.  Furthermore,  we  expect  the  SI  1  responses  at  the  offset  positions  to  become  weaker,  for 
shallow  burial,  due  to  the  cancellation  between  the  diffraction  from  the  top  and  the  bottom  edges.  The 
intensity  of  the  S22  response  (parallel  to  the  plate  orientation)  decreases  in  proportion  to  the 
narrowness  of  the  target’s  top  edge.  Therefore,  for  a  smaller  vertical  plate  this  pattern  will  be  quite 
difficult  to  detect.  Based  on  this  finding,  the  classification  rules  shown  in  Figure  29  need  to  be  modified 
as  shown  in  Figure  40. 

There  were  three  passes  collected  for  the  target  in  cell  HI  1 .  Two  of  the  orthogonal  passes  showed  no 
UXO-like  features.  The  third  pass  showed  features  similar  to  a  vertical  plate  as  shown  in  Figure  44. 

Like  A9,  such  features  led  to  a  UXO-like  classification  at  the  time  of  blind  classification.  On  the  basis 
of  the  modified  classification  rules  below,  HI  1  would  have  been  classified  as  not  UXO-like. 
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Figure  40,  Modified  classification  rules  to  include  features  from  vertical  plates. 

C9  -  New  Classification  Features 

Figure  45  shows  the  classification  features  (ELF  and  ETO)  obtained  from  Item  C9,  a  horseshoe.  What 
is  interesting  about  this  item  is  that  it  is  a  good  resonator  regardless  of  excitation  polarization.  Fiigh 
ELE  was  observed  as  for  an  UXO-Iike  item  without  a  large  depression  angle.  However,  the  ETO  is 
always  the  same  as  the  scan  direction,  which  makes  the  inference  of  target  linearity  suspicious.  Again, 
this  peculiar  feature  was  not  included  in  the  original  classification  features  and  led  to  UXO-like 
classification. 
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Figure  41,  Modified  classification  rules  to  include  features  from  broken  loops  and  vertical 
plates, 

Dl,  D8,  H8  and  H12-  Their  classifications  shown  from 
Figure  46  to 

Figure  49  indicate  horizontal  plate  features  (See  Section  3.2).  The  poor  SCR  of  H12  resulted  in  a  low 
confidence  level. 

H9,  F5,  G9  and  B13  -  The  classification  features  of  these  items  show  high  ELF  regions  but  cannot  be 
confirmed  from  spatial  scattering  pattern  due  to  low  SCR.  Features  of  B13  provide  Figure  50  as  an 
example.  The  other  orthogonal  cut  is  contaminated  by  a  linear  trench  near  the  target  location  and  could 
not  be  used. 
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CIO  &  A6-  Contaminated  by  linear  trenehes.  A  linear  treneh  (near  position  -35  inehes  in  Figure  51  ) 
loeated  near  the  target  makes  the  already  weak  target  response  even  harder  to  diseriminate.  In  faet,  the 
trench  was  mistaken  as  the  target  and  thus  resulting  in  linear  features. 

In  summary,  8  out  of  17  mis-classifications  of  clutter  items  as  UXO-like  could  have  been  avoided 
with  the  inclusion  of  the  additional  classification  rules,  shown  in  Figure  41 .  This  means  a  almost  a 
50%  reduction  of  false-UXO-rate  (FUR)  in  Figure  32. 


fELF  \«  Position  for  E8 
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Estimated  linearity  factor 
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Position  (inches) 

(b)  Parallel  co-polarized  responses. 


S22  E8  Orientation  =  189 


-40  -20  0  20  40 

Position  (inches) 


(c)  Transverse  co-polarized  responses 


Figure  42,  Classification  features  for  Item  E8  (fence  post  cap)  from  survey  line  parallel  to 
target  orientation,  for  fence  post  cap. 
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fELF  vs  Position  for  A9 
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Linearity 

Factor 


(ELF) 


Position  (inches)  Position  (inches) 

(b)  SI  1  Response,  perpendicular  to  plate  orientation  (c)  Transverse  Responses  (S22 

component,  parallel  to  plate  orientation) 


Figure  43  Classification  features  for  Item  A9  (vertical  plate)  obtained  from  a  pass  orthogonal  to 
top  edge  direction  of  the  plate. 
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fELF  vs  Position  for  H1 1 
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(b)  Parallel  eo-polarized  responses. 


S22  H1 1  Orientation  =  90 


-20  0  20  40 

Position  (inches) 

(e)  Transverse  co-polarized  responses. 


Figure  44  Classification  features  for  Item  Hll  (crumpled  box  fin)  from  an  offset  cut. 
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fELF  vs  Position  for  C9  ETO  vs  Position  for  C9  ATV  =  11 


(b)  Estimated  target  orientation  (ETO)  from  parallel  survey  pass. 


fELF  vs  Position  for  C9 


(c)  ELE  from  the  parallel  survey  pass. 


ETO  vs  Position  for  C9  ATV  =  85 


(d)  Estimated  target  orientation  (ETO) 
from  the  orthogonal  survey  pass. 


Figure  45  Classification  features  for  Item  C9  (horse  shoe). 
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ELF  vs  Position  for  D1 
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(b)  Co-polarized  response  parallel  to  pass. 


(c)  Transverse  co-polarized  response. 


Figure  46  Classification  features  for  Item  D1  (M38  box  fin  piece)  from  pass  parallel  to  target 
orientation. 
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fELF  vs  Position  for  D8 
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(e)  Transverse  co-polarized  response. 


Figure  47  Classification  features  for  Item  D8  (M38  box  fin) 
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ELF  vs  Position  for  H8 
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(b)  Co-polarized  response  parallel  to  pass. 


(c)  Transverse  co-polarized  response. 


Figure  48  Classification  features  for  Item  H8  (horizontal  8”  dia,  steel  disk) 
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fELF  vs  Position  for  H1 2 


(a) 

Estimated 

Linearity 

Factor 


(ELF) 


S11  H12  Orientation  =  0 


S22  H12  Orientation  =  0 


Position  (inches) 
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(b)  Co-polarized  response  parallel  to  pass. 


(c)  Transverse  co-polarized  response. 


Figure  49,  Classification  features  for  Item  H12  (box  fin  fragment)  from  pass  parallel  to  target 
orientation. 
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fELF  vs  Position  forB13 


S11  B13  Orientation  =  276 


-20  0  20  40 

Position  (inches) 


(a)  Estimated  Linearity  Factor  (ELF) 


(b)  Co-polarized  response  parallel  to  pass. 


Figure  50  Classification  features  for  Item  B13  (vertical  3”x6”  cyL)  from  survey  pass  parallel 
to  target  orientation 
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(a)  Co-polarized  response  parallel  to  pass. 
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(b)  Transverse  co-polarized  response. 


Figure  51  Classification  features  for  Item  CIO  (vertical  4”x4”  plate)  from  pass  perpendicular  to 
target  orientation. 


69 


5.6  Confidence  Level  and  Signal-to-Clutter  Ratio 


The  confidence  level  was  implicitly  determined  by  signal  intensity,  signal-to-clutter  ratio  and 
consistency  from  multiple-pass  data.  Currently,  the  confidence  levels:  HIGH,  MODERATE  or  EOW 
were  entered  subjectively  by  the  person  performing  the  classification  process.  This  process  started  with 
viewing  pictures  of  features  similar  to  those  shown  in  the  previous  section.  Then  the  operator  went 
through  the  classification  rules  shown  in  Eigure  29  mentally  for  all  passes  available  for  the  target. 

Later,  it  was  found  that  the  SCR  and  signal  intensity  have  very  good  correlation  with  the  confidence 
level  set  subjectively  by  operator’s  impression.  This  means  that  the  confidence  level  can  be  determined 
more  objectively  and  quantitatively  from  the  signal  intensity  and  SCR  in  the  future.  Eigure  52  and 

Eigure  53  plot  the  histograms  of  SCR’s  for  combined 
LOW  &  MODERATE  and  for  HIGH  confidence  levels,  respectively.  These  results  clearly  show  the 
strong  correlation  between  the  confidence  level  and  the  SCR  value.  Note  that  SCR  is  a  function  of 
antenna  and  time  positions.  The  highest  SCR  value  among  all  positions  for  each  item  was  designated  as 
the  SCR  associated  with  that  particular  item.  However,  it  should  also  be  emphasized  that  the 
confidence  level  was  not  determined  by  the  SCR  level  alone.  Whether  features  obtained  from  both 
parallel  and  transverse  pass  showed  consistent  results  (UXO  or  non-UXO)  was  also  one  of  the  factor  in 
specifying  the  confidence  level. 
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Figure  52  Histogram  of  the  SCR  associated 
with  the  low  and  moderate 
confidence  level  cases  (signal  level 
threshold  =  -100  dB  after 
processing). 


Signal  Threshold  =  -lOOdB 
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SCR(dB) 

Figure  53  Histogram  of  the  SCR  associated 

with  the  high  confidence  level  cases 
(signal  level  threshold  =  -100  dB 
after  processing). 
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Chapter  6  Analysis  of  Noise  and  Clutter  Effects  on 
Extracted  Parameters 


6.1  Simulated  Signal  to  Noise  Ratio  (SNR)  Study 

The  study  of  SNR  is  important  for  determining  the  quality  of  the  data  collected,  the  robustness  of  the 
feature  extraction,  and  the  desirability  of  prospective  system  improvements.  In  this  section  we  explore  the 
effects  of  noise  on  the  crucial  parameters  extracted  from  our  GPR  data. 


There  are  many  sources  of  noise  in  a  radar  system,  such  as  noise  from  the  electronics  in  the  network 
analyzer,  thermal  noise  in  the  cable,  and  ambient  noise  in  the  environment.  Noise  is  most  commonly  observed  as 
white  noise,  where  the  amplitude  distribution  at  each  frequency  has  a  Gaussian  shape.  Signal  to  noise  ratio, 

SNR,  is  typically  measured  in  decibels  (dB),  a  logarithmic  scale  value  indicating  the  signal  power  relative  to  that 
of  the  noise 


SNR(dB)  =  10/og 


^signal  ^ 

p  . 

V  noise  J 


(1) 


where  Psignai  and  Pnoise  are  the  respective  powers  of  signal  and  noise.  On  this  scale,  a  factor  of  2  (or  of  one  half) 
change  in  the  power  ratio  implies  a  +(or  -)3  dB  change.  In  general,  note  that  positive  dB  values  imply  a  ratio 
greater  than  one,  while  negative  values  imply  a  ratio  less  than  one.  The  dB  scale  is  also  defined  in  terms  of 
signal  magnitude,  instead  of  signal  power. 


SNR(dB)  =  20 /og 


■^signal 

A 

V  -^noise  J 


(2) 


Because,  generally  speaking,  power  is  proportional  to  amplitude  squared,  these  alternative  definitions  provide 
the  same  results  in  terms  of  dB.  Unless  noted  otherwise,  we  will  use  the  definition  in  (2)  in  what  follows.  For 
reference,  note  that  a  factor  of  10  in  amplitude  ratio  corresponds  to  a  change  of  +20  dB.  Measured  SNR  values 
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quantified  in  this  manner  are  compared  to  desired  values  for  assessing  the  feasibility  of  determining  signal 
content. 


6.1.1  ELF 

In  order  to  understand  how  ELF  and  ETO  [***defmed  above, right?  are  affected  by  noise,  the  following 
idealized  simulation  based  on  an  ideal  linear  (very  elongated)  object  and  a  rotationally  symmetric  object  are 
investigated.  Scattering  from  these  or  other  objects  may  be  expressed  through  an  eigenvalue  matrix,  the  diagonal 
of  which  contains  appropriately  scaled,  intrinsic  scattering  strengths  of  the  object  in  its  principal  orientations,  e.g. 
with  polarization  parallel  to  and  perpendicular  to  its  major  axis.  The  scattering  from  an  infinitely  long  thin  wire 
may  be  represented  in  eigenvalue  matrix.  A,  in  which  one  eigenvalue  is  very  much  dominant,  reflecting  the  fact 
that  scattering  from  a  polarization  perpendicular  to  the  long  axis  is  negligible.  This  implies  a  unitary  ELF. 


Pre 

M 


4 

0  " 

“1 

o' 

0 

0 

0 

^  ELF  =  1 


(3) 


and  post-multiplication  of  A  by  a  rotation  matrix  M  expresses  the  results  under  various  target  orientations, 
is  orthogonal,  i.e.  M  =  M.  For  our  test,  Gaussian  noise  is  added  to  the  scattering  matrix,  S  \ 


Sn  =  S  +  N 


1  0 
0  0 


M  +  N 


where  N  = 


n, 

^3 


n, 

^2 


and  M  = 


cos(y0) 

-sin(y0) 


sin(y0) 

cos(y0) 


(4) 


The  matrix  S  is  the  uncorrupted  signal.  In  the  matrix  N ,  the  components  ni,  n2  and  n3  are  random  numbers  with 

a  Gaussian  distribution,  scaled  by  a  factor  of  10*^  components  nj  and  n2  represent  both  the  internal 

system  noise  and  the  external  ambient  noise.  The  03  component  represents  mainly  the  internal  system  noise  and 
the  noise  produced  from  antenna  element  coupling.  Signal  quantities  are  scaled  relative  to  the  maximum 

eigenvalue,  i.e.  1.  In  the  matrix  M  the  variable  P  is  the  angle  of  the  target  from  a  position  with  axis  parallel  to 

the  signal  polarization.  Thus,  altogether,  the  new  scattering  matrix  Sn  is  the  fundamental  scattering  matrix  S 
with  noise  added.  In  any  given  test  case,  it  is  decomposed  into  its  new  eigenvalues,  whereupon  a  noise-biased 

ELF  is  calculated  from  Sn . 
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Simulations  for  a  thin  wire  were  performed  using  the  procedure  described  above,  with  2000  Monte  Carlo 
realizations.  As  the  SNR  is  varied  from  -10  to  60  dB,  the  added  noise  causes  bias  and  variation  of  the  estimated 
features  (Figure  54  and  Figure  55).  The  rotation  of  the  target,  P,  was  fixed  at  0,  15,  and  45  degrees.  Figure  54 
illustrates  SNR  versus  the  biased  ELF.  As  SNR  declines  to  0  dB  (equal  signal  and  noise  power  and  amplitude), 
the  biased  ELF  approaches  0.45  instead  of  its  true  value  of  1.  The  possible  percent  error  in  the  unbiased  ELF  can 
be  calculated  from  this  graph.  For  example,  the  figure  shows  that  the  SNR  needs  to  be  greater  than 
approximately  20dB  to  produce  a  biased  ELF  greater  than  0.85;  that  is,  to  obtain  a  biased  ELF  with  less  than 
15%  error. 

The  standard  deviation  of  the  biased  ELF  versus  SNR  is  shown  in  Figure  55.  Note  that  the  target 
orientation  has  almost  no  effect  on  the  biased  ELF  for  all  SNR  values.  Flence  accurate  biased  ELF  can  be 
calculated  for  a  target  with  complete  disregard  to  the  orientation  of  the  target.  An  important  revelation  of  this 
exercise  is  that  a  simulated  scattering  matrix  made  up  of  pure  Gaussian  noise  has  an  average  ELF  of 
approximately  0.45.  The  value  0.45  shows  the  natural  offset  or  bias  of  the  ELF  equation.  This  bias  is  due  to  the 
non-linear  operations  used  in  calculating  the  ELF. 


Biased 

ELF 


Figure  54  Biased  ELF  vs,  SNR  Plot  for  P  =  0, 15,  and  45  degrees. 
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Replacing  the  zero  perpendicular  eigenvalue,  X±_,  with  a  one  while  retaining  the  unitary  parallel  eigenvalue 
(above),  simulates  the  biased  ELF  for  a  sphere.  The  eigenvalue  matrix  becomes  symmetric  and  the  ELF 
becomes  zero,  as  it  should  for  a  completely  symmetric  object: 


^// 

0  " 

~\ 

o' 

0 

0 

1 

ELF  =  0 


(5) 
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Figure  55  Standard  Deviation  (STD)  of  ELF  vs,  SNR  for  P  =  0, 15,  and  45  degrees. 


74 


The  eigenvalue  matrix,  A,  is  not  rotated  in  this  test  because  there  is  no  difference  in  the  response  from  different 
target  angles  for  a  completely  symmetric  object.  Gaussian  noise  is  added  to  the  scattering  matrix  in  the  same 

manner  as  for  the  wire  case,  with  the  biased  ELF  subsequently  calculated  from  Sn . 

The  results  of  the  symmetric  object  simulation  are  shown  in  Figure  56  and  Figure  57.  The  former  shows 
that  the  presence  of  noise  causes  the  ELF  to  rise  to  a  value  of  0.53  when  the  SNR  is  around  0  dB.  The  standard 
deviation  of  ELF  is  plotted  versus  the  SNR  in  Figure  57. 
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Figure  56  Biased  ELF  vs,  SNR  for  a  simulated  sphere 


Figure  58  shows  a  comparison  of  the  biased  ELF  for  a  simulated  long  wire  and  a  simulated  sphere.  This  graph 
shows  that  a  target  response  with  SNR  below  approximately  1 0  dB  (the  dark  shaded  region)  would  have  an 
ambiguous  ELF.  This  means  that  the  biasing  of  the  ELF  has  become  so  large  that  it  would  be  very  difficult  to 
determine  a  linear  target  from  a  symmetric  target.  It  appears  that  with  a  SNR  above  approximately  20  dB,  there 
would  be  sufficient  dynamic  range  in  the  ELF  to  determine  the  linearity  of  an  object  reliably.  In  the  "twilight 
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zone"  between  1 0  dB  and  20  dB  SNR,  it  should  often  be  possible  to  determine  the  ELF,  but  possibly  with  some 
ambiguity  or  difficulty. 
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Figure  57  Standard  Deviation  (ELF)  vs.  SNR  for  a  Simulated  Sphere 

The  error  in  the  biased  ELF  for  the  long  wire  and  the  sphere  is  shown  in  Figure  59.  The  error  is 
calculated  by  subtracting  the  theoretical  ELF  value  from  the  simulated  biased  ELF  value.  The  graph  shows  that 
the  error  in  ELF  due  to  noise  is  very  similar  for  both  the  wire  and  the  sphere,  despite  their  gross  differences  in 
geometry.  This  is  quite  useful  for  estimating  the  possible  error  in  ELF  for  all  measured  targets  based  on  their 
measured  SNR.  Next,  we  pursue  a  study  of  the  effects  of  SNR  on  the  ETO. 
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Figure  58  Comparison  of  biased  ELF  vs  SNR  for  a  simulated  long  wire  and  sphere 
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Figure  59  Comparison  of  error  m  ELF  vs  SNR  for  a  simulated  long  wire  and  sphere 


6.1.2  ETO 

The  effect  of  noise  on  the  unbiased  ETO  is  discussed  in  this  section.  The  biased,  i.e.  noise  affected, 

ETO  for  a  wire  is  calculated  from  the  Sn  of  equation  (4)  as  follows.  The  scattering  matrix  Sn  is  decomposed 
into  its  eigenvalues  and  the  associated  eigenvectors  are  determined.  The  direction  in  which  each  eigenvector 
points  is  that  associated  with  its  corresponding  eigenvalue.  Thus,  as  for  the  actual  field  data,  the  ETO  is 
calculated  from  the  arctangent  of  the  dominant  eigenvector  orientation  angles.  Comparison  of  the  biased  ETO  to 
the  actual  target  angle,  P,  reveals  the  biasing  due  to  noise. 

The  biased  ETO  was  calculated  for  various  wire  orientation  angles  and  for  SNR  levels  ranging  from  - 1 0 
to  60  dB.  Figure  60  shows  how  the  unbiased  ETO  feature  is  affected  by  the  presence  of  noise  for  four  target 
angles.  At  a  target  angle  of  zero  degrees  (i.e.  perfect  alignment  with  the  radar  polarization)  the  biased  ETO 
shows  almost  no  dependence  on  SNR,  varying  closely  around  the  correct  value  of  zero.  At  target  angles  of  45, 
60,  and  75  degrees  relative  to  radar  polarization  the  processing  performance  becomes  worse,  the  greater  the 
rotation.  In  reality,  given  the  survey  system  using  at  a  minimum  two  perpendicular  passes  above  a  target 
location,  45  degrees  probably  represents  the  worst  case.  Noting  that  for  P  =  45°  reasonable  ETO  estimates 
appear  for  SNR  less  than  1 0  dB,  we  wee  that  this  parameter  is  less  sensitive  to  noise  than  the  FEE.  The  plot  of 
the  standard  deviation  of  ETO  versus  SNR  is  shown  in  Figure  61  for  the  four  target  angles.  It  illustrates  that  the 
"correct"  values  of  ETO  obtained  at  low  SNR  values  under  zero  rotation  do  not  really  indicate  that  the  system  is 
performing  well  at  that  noise  level.  Rather,  the  system  is  relatively  indiscriminate  there  and  orientations  simply 
tend  to  average  out  to  a  zero  value,  which  happens  to  be  correct. 
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Figure  60  Biased  ETO  vs.  SNR  for  p  =  0,  45,  60,  and  75". 


These  studies  show  how  SNR  can  affect  the  ELF  and  ETO  features.  The  processing  depends 
fundamentally  on  these  features  to  estimate  target  linearity  and  existence  of  a  dominant  orientation,  which  are 


79 


clues  to  UXO-like  geometry.  The  analysis  provides  the  guidance  that  the  measured  data  needs  to  have  an  SNR 
of  at  least  1 0  to  20  dB  to  produce  tolerably  accurate  biased  ELF  and  ETO  values,  with  the  latter  being  less 
sensitive. 

6.2  Measured  Signal  to  Noise  Ratio  (SNR) 

The  analysis  above  has  suggested  that  a  minimum  SNR  level  of  1 0  to  20  dB  is  required  if  we  are  to 
determine  fundamentally  necessary  signal  and  target  features.  The  question  becomes,  how  realistic  is  this  SNR 
threshold?  To  answer  this  question  two  sets  of  data  will  be  examined.  First,  the  method  of  calculating  the  SNR 
from  the  measured  data  will  be  discussed.  Then,  the  distributions  of  SNR  for  the  measured  data  will  be 
presented. 

The  SNR  of  measured  data  is  estimated  in  the  course  of  feature  extraction.  The  magnitude  of  the 
maximum  eigenvalue  in  the  time-domain  is  determined  and  is  then  used  as  a  representative  true  signal  level,  e.g. 
as  zl//  above.  The  noise  level  was  found  by  averaging  a  segment  of  the  negative-time  region  with  the  same 

duration  as  the  signal’s  late-time  region.  The  negative-time  region  was  used  because  it  only  contains  noise.  The 
SNR  was  then  calculated  from  this  representative  signal  level  and  the  noise  level.  For  multi-point  measurements 
the  signal  becomes  a  function  of  position,  with  the  SNR  calculated  at  each  point  in  the  same  way  as  for  the 
single-point  measurements.  The  overall  SNR  value  for  a  site  was  found  by  finding  the  maximum  SNR  for  the 
multi-position  scan. 

Figure  62  shows  the  SNR  values  of  every  target  measured  in  the  previous  Tyndall  AFB  demo, 
superimposed  on  the  average  error  curve  for  the  theoretical  biased  ELF  vs  SNR  (Figure  59).  Note  that  in  this 
case  each  target's  features  were  obtained  using  a  single  spatial  point,  because  at  the  time  of  the  Tyndall  demo 
multi-position  processing  had  not  yet  been  implemented.  For  each  target  the  SNR  from  the  measured  data  was 
used  to  indicate  the  likely  position  on  the  theoretical  error  curve  for  that  target.  Each  error  vs  SNR  value  is 
plotted  on  the  curve  by  a  black  diamond  marker  for  all  153  separate  targets.  Note  that  the  SNR  values  were 
rounded  to  whole  numbers  for  ease  of  plotting;  therefore  there  are  many  overlapping  black  diamond  markers.  In 
any  case,  all  of  the  target  markers  lie  above  the  20  dB  threshold.  In  fact,  all  but  two  targets  have  SNR  values 
greater  than  30  dB,  with  correspondingly  small  associated  error.  This  analysis  strongly  supports  the  inference 
that  noise  is  not  a  limiting  factor  for  this  measurement  set. 
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Figure  62  Tyndall  AFB  Demo  estimated  error  vs  measured  SNR  values  for  each  target 

In  Figure  63  the  SNR  values  for  the  multi-position  measurements  taken  at  Blossom  Point  are  used  to 
locate  the  targets  on  the  same  theoretical  ELF  error  curve  applied  above.  There  were  91  separate  targets.  Again, 
most  of  the  targets  lie  above  the  20  dB  threshold.  Although  at  least  four  targets  show  a  likely  SNR  value  less 
than  20  dB,  once  again  noise  does  not  appear  to  be  a  limiting  factor  for  this  measurement  set.  The  analysis 
allows  us  to  identify  those  targets  with  a  SNR  less  than  20  dB,  which  required  more  careful  examination  and  are 
noted  to  produce  less  confident  classification. 
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Figure  63,  Blossom  Point  estimated  error  vs  measured  SNR  values  for  each  target. 


Figure  64  shows  a  histogram  of  the  SNR  occurrences  for  the  single  point  measurements  at  the  Tyndall 
site.  A  majority  of  the  SNR  values  are  grouped  around  40  dB,  showing  a  very  good  SNR  level  for  the  entire  site. 
Figure  65  shows  a  histogram  of  SNR  occurrences  for  the  multi-position  measurements  of  the  Blossom  Point  site. 
Here  the  SNR  values  are  more  evenly  distributed  between  about  35  to  60  dB.  However,  there  are  8  targets  with 
SNR  levels  less  than  20  dB.  These  8  targets  account  for  less  than  ten  percent  of  the  total  targets.  Again,  it 
appears  that  noise  is  not  a  limiting  factor  for  the  entire  site,  but  could  be  for  a  few  individual  targets,  which  the 
analysis  duly  identifies.  The  other  source  of  interference  is  clutter,  which  will  be  discussed  in  the  next  section. 
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Figure  64.  Distribution  of  occurrences  of  SNR  values  (NSNR)  for  the  Tyndall  Data 
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Figure  65,  Distribution  of  occurrences  of  SNR  values  (NSNR)  for  the  Blossom  Point  Data 


6.3  Simulated  Signal  to  Clutter  (SCR)  Study 

The  study  of  SCR  (dB)  is  also  important  for  determining  the  quality  of  the  data  eollected,  the  robustness 
of  the  feature  extraction,  and  the  desirability  of  various  system  improvements.  In  this  section  we  explore  the 
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effects  of  signal  clutter  on  the  crucial  parameters  extracted  from  our  GPR  data.  Clutter  will  be  defined  as  any 
undesired  signal  that  affects  the  signal  produced  by  the  target  of  interest.  Note  that  we  use  the  term  "signal"  as 
opposed  to  noise.  While  clutter  responses  may  be  complex  and  varied,  they  are  not  random  in  the  sense  that 
system  noise  is,  as  pursued  in  the  previous  section.  Clutter  responses  are  simply  those  reflections  produced  by 
objects  of  no  interest  for  UXO  discrimination,  such  as  rocks,  roots,  cavities,  soil  inhomogeneities,  and 
widespread  small  metallic  scrap.  Individually,  these  signals  are  coherent  in  the  same  way  that  the  desired  signals 
are.  In  most  ground  conditions,  the  average  clutter  levels  in  co-polarization  channels,  Sn  and  S22,  are 
approximately  the  same  and  much  higher  than  that  in  the  cross-polarization  channel,  S21  (by  the  principal  of 
reciprocity,  S21  equals  S12,  so  we  only  define  values  for  one  of  the  cross  polar  channels). 

In  order  to  study  clutter  effects  on  the  extracted  feature  parameters,  simple  simulations  are  performed  in 
which  clutter  responses  are  added  to  the  underlying  target  signal.  Then  the  ELF  and  ETO  for  an  ideal  linear 
object  and  an  ideal  rotationally  symmetric  object  are  obtained  and  evaluated. 


6.3.1  ELF 

In  the  previous  section,  biased  ELF,  ETO,  and  the  scattering  matrix  were  defined.  In  order  to  understand 
how  ELF  is  affected  by  clutter,  one  may  begin  with  a  perfect  linear  scattering  object,  i.e.  an  infinitely  long  thin 

wire.  Then,  as  above,  the  target's  eigenvalue  matrix,  A ,  is  rotated  by  matrix  M  to  simulate  various  target 
orientations.  Clutter  is  added  into  the  scattering  matrix: 


5.  =  m\m  +  c 


where  C  = 


c  c/10 

and  M  = 

cos(y0) 

sin(y0) 

c/10  c 

-  sin(y0) 

cos(y0) 

(6) 


In  the  matrix  M  the  variable  P  is  the  angle  of  the  target  from  an  alignment  parallel  to  the  incident  polarization. 

In  the  matrix  C ,  c  is  a  positive  real  number  scaled  by  a  factor  of  10^  .  Note  that  c  is  simply  a  magnitude, 

not  a  random  variable.  We  do  not  apply  assumed  distribution  to  c,  nor  require  many  Monte  Carlo  realizations  for 

the  analysis.  The  off-diagonal  terms  of  matrix  C  were  also  scaled  by  1/10  to  simulate  a  typical  ratio  between 
the  co-pol  and  cross-pol  channels.  Varying  the  SCR  from  -10  to  60  dB  produced  different  values  of  c.  The 
signal  eigenvalue  matrix  was  either  that  for  the  long  wire  or  sphere: 
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After  the  clutter  signal  was  added  to  the  system,  the  new  scattering  matrix,  Sc ,  was  decomposed  into  it’s 
eigenvalues  and  the  ELF  determined. 


Figure  66  and  Figure  67  show  that  the  added  clutter  causes  a  bias  in  the  estimated  features.  For  the  wire  ( 
Figure  66),  equal  signal  and  clutter  magnitudes  (0  dB  SCR)  produce  a  biased  ELF  of  only  0.33,  relative  to 
the  "true"  value  of  one.  Overall,  the  SCR  must  reach  about  20  dB  before  the  biased  ELF  error  is  only  20%  or 
less.  Interestingly,  there  was  almost  no  discernable  difference  in  the  biased  ELF  curves  for  various  target  angles, 
P,  therefore  only  the  zero  degree  case  is  shown.  This  supports  the  position  that  accurate  ELF  can  be  calculated 
for  a  target  with  complete  disregard  to  the  orientation  of  the  target,  as  observed  above  in  the  SNR  analysis. 

For  the  sphere,  the  ideal,  noise  and  clutter-free  ELF  is  zero.  Figure  67  shows  that  the  biased  ELF  rises  to  a 
value  of  0.05  when  the  SCR  declines  to  0  dB.  This  shows  that  the  ELF  estimation  for  a  symmetric  object  is  not 
effected  by  the  presence  of  symmetric  clutter.  Figure  68  shows  the  error  in  biased  ELF  for  both  the  sphere  and 
the  long  wire.  The  two  error  curves  will  be  used  below  in  discussion  of  measured  targets. 
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Figure  66.  Biased  ELF  versus  SCR  for  a  simulated  long  wire 
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Figure  67.  Biased  ELF  versus  SCR  for  simulated  sphere 
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Figure  68,  Comparison  of  error  in  ELF  for  simulated  long  wire  and  sphere. 


6.3.2  ETO 

The  effect  of  simulated  clutter  on  the  unbiased  ETO  is  discussed  in  this  section.  As  before,  the  biased 

ETO  is  calculated  from  Sc  in  equation  (6).  The  scattering  matrix  Sc  is  decomposed  into  its  eigenvalues  and 
eigenvectors;  then  the  ETO  is  calculated  as  the  arctangent  of  the  eigenvector  angles.  The  difference  between  the 
biased  ETO  and  the  target  angle,  P,  is  the  biasing  due  to  clutter. 
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The  biased  ETO  was  calculated  for  a  long  wire  with  various  rotation  angles  relative  to  alignment  with 
the  transmitted  polarization.  The  SCR  level  ranged  between -10  to  60  dB.  Figure  69  shows  how  the  unbiased 
ETO  feature  is  affected  by  the  presence  of  clutter  for  four  target  angles,  by  plotting  the  error  in  biased  ETO 
versus  SCR.  At  a  target  angle  of  45  degrees  the  ETO  shows  no  dependence  on  SCR.  At  target  angles  of  0 
(upper  curve),  75,  and  90  degrees  the  ETO  error  does  not  become  entirely  negligible  until  the  SCR  is 
approximately  10  dB.  However,  the  ETO  error  is  tolerable  at  an  SCR  level  of  about  0,  and  is  only  15%  for  an 
SCR  at  the  worst  value  of  -10  dB.  That  is,  the  ETO  error  is  in  general  not  significant. 


ETO 

Error 

(deg) 


-10  0  10  20  30  40  50  60 

SCR(dB) 


Figure  69  ETO  error  in  degrees  vs  SCR, 

These  calculations  support  the  view  that  the  measured  data  needs  to  have  an  SCR  greater  than  about  20  dB 
in  order  to  produce  a  valid  estimate  of  the  EEF.  The  ETO  is  much  less  sensitive,  with  very  small  clutter¬ 
generated  error  overall.  Thus,  as  above  in  the  examination  of  SNR  effects,  the  signal  contamination  is  most 
likely  to  limit  determination  of  EEF  but  not  ETO. 
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6.4  Measured  Signal  to  Clutter  Ratio  (SCR) 


The  signal  to  clutter  ratio  (SCR)  for  the  measured  data  from  the  Tyndall  and  Blossom  Point  sites  will  be 
presented  in  this  section.  First,  the  method  of  calculating  the  SCR  from  the  measured  data  will  be  discussed. 
Then,  the  distributions  of  SCR  for  the  measured  data  will  be  presented. 

The  SCR  of  the  measured  data  was  calculated  after  the  feature  extraction.  The  magnitude  of  maximum 
eigenvalue  in  the  time-domain  was  found  from  the  measured  data.  The  eigenvalue  was  used  for  the  signal  level. 
The  clutter  level  was  found  by  using  the  response  of  an  empty  site  and  extracting  the  eigenvalues  from  its  late¬ 
time  region.  The  empty  site  late-time  region  was  set  equal  to  the  late-time  region  of  each  target.  This  ensured 
corresponding  late-time  regions,  with  similar  depths.  The  maximum  empty  site  eigenvalue  was  set  as  the  clutter 
level.  The  SCR  was  then  calculated  from  the  signal  level  and  the  clutter  level  of  each  site. 

Figure  70  shows  a  single  position  time-domain  response  for  a  measured  target.  The  highlighted  region  is 
the  selected  late-time  region  (12-26  ns)  used  for  feature  extraction  and  the  calculation  of  the  signal  level. 
Figure  71  shows  a  single  position  time-domain  response  for  a  measured  empty  site.  The  highlighted  region 
corresponds  to  the  late-time  region  (12-26  ns)  selected  for  the  target  in  Figure  70  and  is  used  to  calculate  the 
clutter  level. 
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Figure  70  Target  response  S  and  selected  late  time  region  (signal) 
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Figure  71  Empty  site  response  and  corresponding  selected  time  region  (clutter) 


For  multi-point  measurements  the  SCR  beeame  a  funetion  of  position.  The  SCR  is  ealeulated  the  same  as 
in  the  single-point  measurements  for  eaeh  antenna  position.  The  overall  SCR  value  for  a  site  was  found  by 
finding  the  position  with  the  maximum  elutter  signal  level. 

Figure  72  shows  the  "waterfall"  plot  of  a  target  measured  using  the  multi-position  survey  method  at 
Blossom  Point.  The  two  white  lines  superimposed  on  the  plot  represent  the  seleeted  late-time  region  for  the 
given  target.  This  region  was  used  to  ealeulate  the  signal  level  for  this  target.  The  same  late-time  region  lines 
are  superimposed  on  Figure  73,  the  waterfall  plot  of  an  empty  BP  site.  This  region  was  used  to  ealeulate  the 
elutter  level.  The  empty  site  has  the  same  frequeney  filtering  and  smoothing  windows  applied  to  it  as  the  target 
response  had.  This  is  done  to  try  to  ealeulate  a  fair  representation  of  the  elutter  level. 
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Figure  72  Target  (MK-23)  site  (B2)  waterfall  plot  (signal) 
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Figure  73  Empty  site  waterfall  plot  (clutter) 

In  Figure  74  the  SCR  values  of  every  Tyndall  demo  target,  measured  using  the  single  point  method,  are 
superimposed  on  the  two  theoretical  biased  ELF  error  discussed  above.  That  is,  each  SCR  value  is  plotted  on  the 
ELF  curve  with  a  black  diamond  marker  on  the  long  wire  error  and  a  black  fdled  circle  for  the  sphere  error. 

There  were  1 53  separate  targets  measured  at  Tyndall  AFB  (Appendix  A).  It  should  be  noted  that  the  SCR  values 
were  rounded  to  the  nearest  one  tenth  for  ease  of  plotting;  therefore  there  are  many  overlapping  black  diamond 
and  circle  markers.  One  can  see  that  all  of  the  targets  lie  between  5  and  25  dB.  A  majority  of  the  targets  lie 
below  the  20  dB  threshold  discussed  earlier.  This  shows  that  there  will  possibly  be  errors  in  the  ELF  and  ETO 
parameters  due  to  clutter. 
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Figure  74  Theoretical  error  vs  SCR  for  each  target  in  the  Tyndall  data 

A  histogram  of  the  measured  SCR  for  the  single  point  measurements  at  the  Tyndall  site  is  shown  in 
Figure  75.  It  can  be  seen  that  a  majority  of  the  measured  SCR  values  are  grouped  around  10  dB,  showing  a 
moderately  poor  SCR  level  for  the  entire  site.  It  appears  that  clutter  is  a  limiting  factor  for  this  entire  site.  The 
question  becomes,  “How  do  we  reduce  clutter?”  One  possible  answer  can  be  seen  in  the  multi-position 
measurements  taken  at  Blossom  Point. 
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Figure  75.  Number  of  SCR  values  (NSCR)  vs  SCR  for  Tyndall  data 


In  Figure  76  the  measured  SCR  values  for  the  multi-position  measurements  taken  at  Blossom  Point 
Proving  Grounds  are  superimposed  on  the  theoretical  biased  ELF  plot.  These  SCR  values  were  calculated  using 
the  same  procedure  used  to  calculate  the  single  point  SCR  from  the  Tyndall  site.  The  SCR  levels  for  the  targets 
range  from  -10  to  20  dB.  It  appears  that  all  of  the  targets  lie  below  the  20  dB  threshold  that  was  discussed 
previously.  This  allows  an  apples-to-apples  comparison  with  the  Tyndall  data,  strongly  suggesting  that  there  is 
an  overall  lower  SCR  level  for  the  Blossom  Point  site.  A  histogram  of  the  SCR  for  the  Blossom  Point  site  is 
shown  in  Figure  77,  with  SCR  values  clustered  between  mainly  0  to  10  dB.  It  appears  that  clutter  is  limiting 
factor  for  the  entire  site. 


93 


1 

0.9 
0.8 
0.7 

Theo-  j 
retical  ^  0.6 

D 

error  o 

5  0.5 

D 

:  0.4 

D 

j  0.3 
0.2 
0.1 
0 

-10  0  10  20  30  40  50  60 

SCR(dB) 

Figure  76  Theoretical  error  in  Blossom  Point  Data  vs  SCR  values  for  each  target;  SCR 
calculated  using  same  method  as  at  Tyndall  for  each  antenna  position. 
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Figure  77  Number  of  SCR  level  occurrences  (NSCR)  for  Blossom  Point  Data,  vs  SCR,  calculated 
using  same  method  as  at  Tyndall,  for  each  antenna  position. 

Now  let’s  look  at  the  measured  SCR  values,  when  calculated  using  the  multi-position  processing  method, 
for  the  multi-position  measurements  taken  at  Blossom  Point.  Values  are  superimposed  on  the  theoretical  biased 
ELF  plot  in  Figure  78.  One  can  see  that  the  SCR  levels  for  the  targets  range  from  5  to  60  dB.  It  appears  that  a 
majority  of  the  targets  lie  above  the  20  dB  threshold  that  was  discussed  previously.  Furthermore,  this  plot  shows 
that  there  are  several  targets  with  a  SCR  value  less  than  20  dB.  A  histogram  of  the  SCR  for  the  multi-position 
measurement  Blossom  Point  site  is  shown  in  Figure  79.  Relative  to  the  distribution  in  Figure  77  the  SCR  values 
are  much  more  evenly  distributed,  now  between  about  15  to  55dB.  However,  there  are  approximately  19  targets 
with  a  SCR  levels  less  than  20  dB.  These  19  targets  account  for  20  percent  of  the  total  targets.  It  appears  that 
clutter  is  not  a  limiting  factor  for  the  entire  site,  but  it  is  still  a  concern.  This  demonstrates,  in  part,  one  of  the 
advantages  of  multi-position  measurements  over  the  single-point  measurements.  The  main  difference  was  the 
background  subtraction.  Single  point  background  subtraction  is  performed  by  selecting  one  position  as  target 
location  and  another  as  background.  While  multi-position  background  subtraction  uses  every  position  as  the 
target  and  uses  an  ensemble  average  of  all  the  positions  as  a  background. 
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Figure  78  Theoretical  error  in  Blossom  Point  data  vs  SCR  values  for  each  target. 
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NSCR 


Figure  79  Number  of  occurrence  of  SCR  levels  (NSCR)  for  multi-position  treatment  of  Blossom 
Point  Data. 
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Chapter  7  Conclusions 


The  old  lanes  at  the  Blossom  Point  site  eontain  many  vertieal  UXO-like  items.  This  ordnanee 
orientation  is  in  general  very  diffieult  for  GPR  operations.  First,  a  radar  direetly  above  the  target  exeites 
only  weak  seattering.  Although  the  seattering  ean  be  enhaneed  by  offsetting  the  radar  to  provide 
oblique  subsurfaee  ineidenee,  a  longer  propagation  path  results  in  more  soil  absorption.  During  the 
Demo  I  in  this  projeet  (Tyndall  AFB),  it  was  already  determined  that  the  elassifieation  performanee  for 
UXO-like  items  with  large  inelination  angles  was  poor.  Therefore,  the  faet  that  only  one  vertieal 
ordnanee  item  was  missed  during  the  blind  test  at  Blossom  Point  is  indeed  a  remarkable  improvement. 
Sueh  improvement  was  aehieved  for  two  reasons.  First,  the  eomprehensiveness  of  the  target  set  at 
Tyndall  provided  invaluable  lessons  as  far  as  finding  out  the  sourees  of  missed  UXO  and  false 
elassifieation.  Seeondly,  based  on  the  Tyndall’s  findings  and  the  radar  knowledge  of  our  team, 
effeetive  improvements  in  system  eonfiguration,  measurement  approaeh,  feature  extraetion  algorithms 
and  elassifieation  eriteria  were  made  before  the  Demo  II  at  Blossom  Point  Demo.  In  part  the  new 
approaeh  involved  making  multiple  passes,  in  different  direetions,  over  essentially  all  targets.  Single 
passes  over  some  targets  were  assoeiated  with  poorer  performanee  in  the  previous  demo. 

The  new  multiple-position,  multiple  pass  approaeh  proved  to  be  useful  in  diseriminating  UXO-like 
items  from  elutter  items  for  ordnanee  with  large  inelination  angles,  by  providing  spatial  variation  of 
target  features  and  temporal-spatial  seattering  patterns.  Canonieal  seattering-pattem  features  for 
eondueting  spheres,  plates  and  UXO-like  items  were  defined  and  applied  for  elassifieation.  Some 
unique  site-speeifie  elutter  items  sueh  as  horseshoes  and  barbed  wire  bundles  also  showed  interesting 
linearity  and  seattering  features  that  are  different  from  a  UXO-like  item,  elassifieation  rules  that  utilize 
estimated  linear  faetor  (ELF),  estimated  orientation  (ETO)  and  seattering  patterns  were  developed  and 
refined  based  on  these  eanonieal  objeets.  Currently,  the  elassifieation  proeedure  was  performed  by 
training  an  operator  to  reeognize  visuals  in  exeeuting  the  elassifieation  rules.  While  requiring  “human 
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in  the  loop”  processing,  the  training  required  is  not  formidable.  Automatic  classification  can  be 
implemented  with  further  development. 

The  classification  performance  on  the  new  lanes  at  Blossom  Point  was  much  poorer  than  on  the  old 
lanes.  We  conclude  that  this  was  related  to  three  facts:  (1)  Many  small  items  were  buried  at  much 
greater  depth  than  the  similar  items  buried  in  the  old  lanes;  (2)  The  trenches  caused  by  the  burial  of  the 
targets  were  more  prominent  than  on  the  old  lanes.  There  were  also  visible  linear  trenches  existing  in 
the  whole  Blossom  Point  site,  perhaps  caused  by  previous  burial  or  extraction  of  drainage  or  irrigation 
system;  (3)  more  inclination  variations  in  the  new  lanes.  Overall,  one  notes  that  UXO-like  items  with  a 
large  depression  angle  (>  45°)  and  depth  greater  than  0.75  meters  showed  poor  signal  to  clutter  ratio 
and  were  much  more  difficult  to  classify,  under  the  soil  conditions  at  the  site.  The  combination  of  weak 
scattering  and  greater  soil  absorption  resulted  in  weaker  radar  response.  The  sensitivity  of  the  radar 
itself  is  greater  than  90  dB  and  is  capable  of  detecting  such  weak  responses.  However  the  clutter  level 
(approximately  70dB  below  the  input  power)  resulting  from  ground  scattering  (surface  &  subsurface)  is 
a  major  limitation.  Increasing  transmitted  power  will  only  increase  both  the  signal  and  clutter  levels 
and  does  not  improve  the  SCR.  In  the  very  least,  it  may  be  advisable  in  some  circumstances  to 
segregate  signal  regions  for  processing  so  that  near-surface  and  deeper  reflections  are  treated 
separately,  each  ear-marked  so  that  confidence  can  be  assigned  accordingly. 

Some  possible  approaches  for  improving  the  SCR  include:  (1)  bistatic  measurement  configuration 
that  uses  separate  transmitting  and  receiving  antennas;  (2)  further  clutter  reduction  using  signal 
processing  methods;  (3)  smart  radiation  pattern  control  of  the  antenna;  (4)  smoothing  the  ground 
surface.  The  effectiveness  of  the  Approach  (1)  depends  on  whether  the  ground  clutter  scatters  more 
weakly  into  the  bistatic  direction  compared  to  the  backscattering  direction.  This  is  yet  to  be  evaluated, 
and  may  depend  on  the  type  of  clutter.  Similar  to  the  adaptive  smoothing  algorithm  that  was  developed 
here  to  greatly  improve  the  SCR  utilizing  the  different  propagation  delays  between  the  desired  target 
and  undesired  clutter  sources,  other  processing  algorithms  can  be  developed  to  utilize  other 
distinguishing  features  between  the  target  and  different  types  of  clutter.  A  majority  of  the  clutter  is 
caused  by  ground  surface,  not  only  from  the  spot  underneath  the  antenna  but  also  from  spots  near  the 
antenna.  Current  antenna  design  has  a  relatively  broad  radiation  pattern  in  order  to  have  effective 
illumination  at  the  target  from  offset  positions.  This  design  also  makes  it  vulnerable  to  the  scattering 
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from  the  ground  surrounding  the  antenna.  A  more  desirable  design  is  to  have  a  narrow  radiation  pattern 
but  with  controllable  beam  direction  as  illustrated  in  Figure  80.  Although  this  is  a  very  sound  design 
conceptually,  it  is  a  big  challenge  for  GPR  antenna  design  due  to  the  frequency  range  and  limited 
antenna  size.  A  related  alternative  is  to  use  (perhaps  in  combination  with  the  system  already  deployed) 
a  special  purpose  antenna,  with  fixed  beam  orientation  but  narrower  beamwidth  and  oblique  subsurface 
direction.  Development  of  an  antenna  with  such  a  design  has  been  proposed  for  a  SERDP  project. 
Approach  (4)  might  seem  unrealistic,  however  at  present  it  may  be  the  most  effective  way  of  reducing 
the  surface  clutter. 


Air 

\\\  u// 

Soil 


Figure  80  Conceptual  smart  antenna  design  with  narrow  beam  patterns  (red)  and  controllable 
beam  direction,  shown  oriented  in  selected  directions. 

The  Blossom  Point  Site  data  indicated  that  signals  collected  from  a  transverse  (second)  pass 
centered  at  target  location  were  not  useful  for  ordnance  with  a  significant  inclination.  A  more  useful 
transverse  pass  would  be  collected  from  the  offset  position.  Figure  81  and  Figure  82  give  two  examples 
of  buried  UXO’s  that  have  depression  angles  of  30  degrees.  The  data  were  collected  along  a  UXO 
direction  centered  at  target  locations  with  the  antenna  polarized  parallel  to  the  UXO’s  azimuthal  angle 
(horizontal  orientation).  Both  data  clearly  show  the  depression  angle  from  the  tilted  reflection 
responses  that  are  offset  from  the  center,  associated  with  strong  backscatter  along  the  incident  angles 
shown  in  Figure  83.  The  previous  approach  collected  the  transverse-pass  (out  of  paper)  data  centered  at 
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the  target,  i.e.  through  position  0.  As  one  can  see,  the  responses  observed  at  the  center  position  are 
actually  quite  weak.  If  one  selected  the  center  position  of  the  transverse  pass  based  on  the  parallel-pass 
data,  e.g.  the  position  at  30  inches  in  Figure  81,  then  a  better  responses  can  be  obtained  from  the 
transverse  pass  with  a  likely  improvement  in  classification  capability.  This  new  approach  is  being 
tested  at  the  OSU/ESL  test  field. 

Targets  contained  in  Tyndall  and  Blossom  Point  sites  represent  only  a  small  portion  of  real-world 
items.  Therefore,  conducting  more  measurements  at  different  sites  containing  different  ground 
conditions,  UXO  types,  and  clutter  items  is  very  important.  This  is  underscored  by  the  lessons  learned 
and  improvements  in  survey  approach  and  processing  that  were  developed  in  response  to  each  of  these 
Demos. 


Position  (inches) 


Figure  81  GPR  data  collected  parallel  to  the  UXO  (2,75”,  Position  19)  orientation  that  has  a  30- 
degree  depression  angle. 
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113(0°, 30°) 


Position  (inches) 


Figure  82  GPR  data  collected  parallel  to  the  UXO  (81mm,  Position  113)  orientation  that  has  a  30- 
degree  depression  angle. 
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Figure  83  Stronger  scattering  for  a  tilted  UXO  from  offset  position. 
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Appendix  A  Archiving 

All  the  raw  and  processed  data  will  be  submitted  in  a  self-documented  CD-ROM.  A  couple  of 
MATLAB  programs  that  allow  the  user  to  read  the  data  (time  and  frequency)  and  view  the 
classification  features  will  also  be  provided. 

A.l  GPR  Data  Files  -  *.cdt 

The  data  for  system-calibrated  frequency-domain  radar  data  was  stored  in  Binary-format  files 
called  “aydddaa.cdt”,  where  “a”  is  from  A  to  Z  for  file  ordering.  The  letter  “y”  indicates  the 
last  digit  of  the  year.  For  example,  “0”  represents  the  year  of  2000.  The  three-digit  number, 
“ddd”  ,  indicates  the  Julian  date  when  the  data  was  stored.  Each  file  contains  a  two-dimensional 
data  taken  at  multiple  target  locations.  The  first  frequency  (10  MHz)  data  was  stored  in  the  first 
row,  the  second  frequency  (12  MHz)  data  was  stored  in  the  second  row,  . . .,  etc.  Table  9 
illustrates  the  new  format  of  multiple  GPR  data.  Data  collected  at  each  position  is  contained  in 
one  vertical  block  that  is  compose  of  six  columns  corresponding  to  the  real  and  imaginary  parts 
of  Si  1,  S21  and  S22,  respectively,  as  illustrated  in  Table  10.  For  each  radar  file,  ‘*.cdt”,  there  is 
an  associated  comment  file  called  “*.txt”  to  store  the  system  information  and  comments.  There 
is  also  an  associated  processed  file  called  “*.mat”  file  that  stores  the  processed  results  and  all 
processing  parameters.  All  of  these  files  will  be  available  in  a  CD-ROM  after  this  submission  of 
this  report. 

Table  9  Format  for  multiple-position  GPR  data. 


Data(fi,Xi) 

Data(fi,X2) 

Data(fi,Xn) 

Data(f2,Xi) 

Data(f2,X2) 

Data(f2,Xn) 

Data(fm,Xi) 

Data(fm,X2) 

D  3.t3.(  fni,Xn) 

An  example  of  data  format  for  Data(fm,Xn)  is  shown  below. 
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Table  10  Data  format  for  a  single  frequency  at  a  single  position. 

R6(Sii)  lm(Sii)  R6(S2i)  lm(S2i)  R6(S22)  1^(822) 


"ReQ"  and  "Im()"  indicate  the  real  and  imaginary  parts,  respeetively  (eombined,  these  provide 
the  amplitude  and  phase). 

A.2  Comment  Text  Files  -  *.txt 

The  eomment  text  files  eontain  information  about  measurement  eonditions  (i.e.  position, 
direetion,  etc.)  and  any  comments  the  user  entered  during  the  measurements.  For  example,  the 
eomment  file  a0348ag.txt  is  printed  below: 

Start  Time:  511/Stop  Time:  519/Target  #:  A2/Start  f:  10/Stop  f:  410/Num  of  Points:  201/Antenna 
Orientation:  92/Num  of  Scans:  41/Antenna  Position:  xstart:  -60/deita  x:  3/yoffset:  0/Reiative 
Permitivity:  20/User  Comments:  /$ 

This  is  a  typieal  eomment  file.  It  eontains  the  time,  in  minutes  from  midnight,  at  whieh  the 
sean  was  started  and  finished.  It  also  eontains  the  information  neeessary  to  reereate  the 
frequeney  veetor  and  the  position  veetor,  i.e.  number  of  positions  of  data  taken  and  the  spaeing 
between  them. 

A.3  Processed  Files  -*.mat  Files 

The  proeessed  data  is  saved  in  *.mat  format  (a  Matlab  file  format),  in  whieh  the  following 
variables  are  saved. 

Variables  of  Processed  Data  Results 

ELF  -  Estimated  Linear  Factor  vs.  position 

fELF  -  Frequency  Estimated  Linear  Factor  vs.  position 

melf  -  Mean  of  ELF  and  fELF 

ETO  -  Estimated  Target  Orientation  vs .  position 

CNR  -  Complex  Natural  Resonance  and  Damping  vs.  position 

ETL  -  Estimated  Target  Length  vs.  position 

Yla  -  Position  vector 

SNR  -  Estimated  Signal  to  Noise  Ratio  vs.  position 

SCR  -  Estimated  Signal  to  Clutter  Ratio  vs.  position 

ATV  -  Antenna  Orientation 
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Variables  of  How  Data  was  Processed 

SELECTION  -Data  Channel  used  for  CNR  extraction 
ftype  -  Type  of  frequency  domain  filter 

ucfreq2  -  User  Centered  frequency  for  the  adaptive  bandpass  filter 

npoint  -  Number  of  points  in  the  adaptive  bandpass  filter 

nx  -  Position  vector  for  late-time  region  selection 

ntmax  -  First  Time  position  for  late-time  region  vs .  x-position 

nuitiTdiff  -  Number  of  points  in  late-time  region  vs.  x-position 

TO  -  Variable  for  slope  gain 

TM  -  Variable  for  slope  gain 

gainfac  -  Variable  for  slope  gain 

imgT  -  Variable  for  adaptive  smoothing 

imgindT  -  Variable  for  adaptive  smoothing 

imgindx  -  Variable  for  adaptive  smoothing 

imgN  -  Variable  for  adaptive  smoothing 

SIDE  -  Variable  for  adaptive  smoothing 
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Appendix  B  Electrical  Properties  of  the  Blossom 


Point  Soil 


The  soil  electrical  properties,  permittivity  and  conductivity  at  40  MHz  and  60  MHz,  were 
measured  at  different  depths  using  the  OSU/ESL  soil  probe  as  shown  in  Figure  84 


Figure  84  OSU/ESL  soil  probe  for  permittivity  and  conductivity  measurement  at  40  MHz 
and  60  MHz, 
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24  feet  east  of  El  2  center 


Figure  85  Soil  electrical  property  measured  on  12/10/2000  at  40  MHz, 


24  feet  E  of  E1 2  Center 


Figure  86  Soil  electrical  property  measured  on  12/10/2000  at  60  MHz, 
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30  feet  E.  of  E12 


Figure  87  Soil  electrical  property  measured  on  12/11/2000  at  40  MHz, 


30  feet  E.  of  El  2 


Figure  88  Soil  electrical  property  measured  on  12/11/2000  at  60  MHz, 
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Middle  of  A8  &  18 


Figure  89  Soil  electrical  property  measured  on  12/12/2000  at  40  MHz, 


Middle  of  A8  &  18 


Figure  90  Soil  electrical  property  measured  on  12/12/2000  at  60  MHz. 
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Appendix  C  BLOSSOM  PT  DEMO  LOG  Excerpt 


BLOSSOM  POINT,  12/13-14/2000,  SET  1  SCHONSTEDT  INEO 


North 


t 


15  B  -  dipole  N-S 

15  C  -  ~  shallow,  dip  on 

15  D  -  shallow  dip 

15  E  -  very 
shallow 

15  A  -  the  reference 

re-examination,  S 

~  NW  -  SE 

14  A  -  shallow  N-S 
dip 

14  B  -  shallow 

14  C  -  shallow,  N-S 
dipole 

14  D  -  shallow, 
symm 

14  E  -  broad, 
south 

pc-d 

fmissed  shallow?! 

13  A  -  strong 

13  B  -  shallow 

13D-C 

13E-X 

12  A  -  very  shallow 

12  B  -  dipole,  Nal 

X 

12  D  -  NR 

12  E  -shallow  dip 
slight  indie  of  Nal 

11 A -strong  shallow 
dip,  aligned  with 
marked  directions 

11  B  -  shallow, 
strong 
dip  ~  N-S 

11 C  -  shallow,  dip, 
aligned 

11  D  -  weak,  S 

llE-X 

X 

10  B  -  shallow 

N-S  dipole 

10  C  -  very  shallow 

lOD-X 

10  E  -  S/SE 

9  A  -  very  shallow 

9  B  -  shallow 

9  C  -  shallow,  dip 
~  aligned 

9D-C 

9  E  -  shallow 

8  A  -  very  shallow 

X 

8  C  -  shallow,  aligned 
dipole 

8  D  -  strong, 
sombrero? 

8  E  -  shallow, 
slight 

indie  of  dip 

7A-C 

7  B  -  shallow 

X 

7D-C 

7  E  -  centered 

6A-NR 

X 

6  C  -  shallow,  some 
dipole 

X 

6E-X 

5  A  -  fairly  shallow,  S 

5  B  -  a  bit  shallow, 

S 

X 

5  D  -  shallow 

5E-X 

4  A  -  South 

4B-C 

4  C  -  shallow 

4  D  -  broad 

4E-X 

3A-C 

3B-C 

3C-NR 

3  D  -  broad,  south 

3E-NR 

2  A- strong  dipole  in 
all  directions 

2  B  -  shallow. 

South, 
slight  dipole 

2  C  -  dipole,  N-S 

2  D  -  shallow, 
sombrero 

2  E  -  broad,  south 

1 A  -  shallow 

X 

1 C  -  weird, 
heterogeneous,  S 

1 D  -  shallow,  dip 

1 E  -  shallow,  dip 

C  =  centered,  azimuthally  symmetric,  positive;  dip  =  dipole  (sign  change) 

S,  NE  etc  =  south,  northeast  etc  offset  NR  =  no  significant  signal 

Nal  =  not  aligned  with  marked  principal  directions  X  =  no  target 

highlight  =  calibration  targets;  NOTE:  The  flag  for  13  A  was  put  on  12  A  by  mistake  for  at  least  part  of 


the  time. 


no 


BLOSSOM  POINT,  12/13/2000,  SET  4  SCHONSTEDT  INFO 


North 


15 1  -  C,  broad 

14 1  -  South 

13 1  -  a  bit  shallow. 
South 

12  H  -  SE 

12 1  -  C,  broad 

IIH-C 

11 1  -  C,  broad 

10  H  -  broad 

10 1  -  C,  broad 

9  G  -  faint 

9H-SW 

91-C 

8G-NR 

8  H  -  South 

81 -NR 

7F-C 

7G-C 

7H-C 

71 -NR 

6F-C 

6  G  -  broad,  faint 

6  H  -  South 

61 -NR 

5F-NR 

5  G  -  moderately 
broad 

5  H  -  shallow,  NE 

51 -NR 

4G-NR 

4  H  -  South 

41-C 

C  =  centered,  azimuthally  symmetric,  positive; 
dip  =  dipole  pattern  (sign  change) 

S,  NE  etc  =  signal  center  offset  south,  northeast  etc 
NR  =  no  significant  signal 

Nal  =  not  aligned  with  marked  principal  directions 
blank  cell  =  no  target 
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Appendix  B:  Ground  Truth  Data 


The  following  information  about  the  ground  truth  data  for  the  Blossom  Point  site  was  provided 
by  NRL  after  the  submission  of  the  processing  results  to  ESTCP.  Figure  91  shows  the  location 
designation  of  the  targets  buried  at  the  Blossom  Point  site.  Only  the  targets  buried  in  lanes  from 
A  thru  I  were  measured  by  CRREL/OSU  team.  The  new  GPR  targets  were  buried  in  lanes  from 
F  to  I.  The  ground  truth  of  these  targets  are  shown  by  the  target  ID  shown  in  Table  1 1  and 
Figure  92. 
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Figure  91  Location  designation  of  the  Blossom  Point  site.  Only  lanes  A~I  were  measured  by 
OSU/ESL  GPR. 
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Table  11  The  target  ID  (Ground  Truth)  of  the  targets  buried  in  Lanes  A~E, 


15 

GO 10250 

LO 10250 

BO 10200 

AO  10050 

KOI 0080 

14 

HO 10052 

MO  1025 

JO 10252 

C010350 

E010750 

13 

P010750 

DOlOlOO 

E010500 

N010350 

12 

B020100 

L020150 

Q070100 

K020050 

11 

Q040030 

D040100 

Q090100 

E021000 

10 

M02025 

H020050 

E030154 

9 

J020250 

E020350 

Q120050 

D020250 

Q060050 

8 

A020100 

K030054 

Q100150 

Q150080 

7 

C030500 

F030750 

E040250 

D030250 

6 

0030030 

B030100 

5 

C030350 

Q140230 

M03015 

4 

N040500 

E030350 

B040200 

E041000 

3 

C040500 

N030500 

Q010050 

E050254 

Q020050 

2 

Q080230 

M04015 

QllOlOO 

E060150 

P020750 

1 

N040350 

E040500 

Q050100 

Q130150 

A 

B 

C 

D 

E 

Item  Designations  are  of  the  form  AxxDDDAI  where  A  =  Item 
Designator,  xx  =  Serial  Number,  DDD  =  depth  in  em,  A  =  azimuth  in 
1/16  of  a  cirele  (0  =  0  deg  ,  2  =  45  deg,  etc.),  and  I  =  inclination  in  the 
same  units  (0  =  parallel  to  the  surface,  4  =  90  deg  nose  down,  etc.) 
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Figure  92  Target  ID  of  the  new  GPR  targets  buried  in  lanes  F~I. 


Table  12  Description  of  the  targets  buried  in  lanes  A~E. 


I 

Descript, 

I 

Descript. 

I 

Descript. 

I 

Descript 

I 

Descript, 

A 

B 

C 

1i/2"x12'' 

steel 

D 

3"x6" 

steel 

E 

3"xl2'' 

steel 

F 

3"x24" 

steel 

G 

161b 

shotput 

H 

4"x4"xy4" 
steel  plate 

J 

8"x8"xy4 
"  steel 

K 

ii/2"x6"xy4 
"  steel 

L 

3"xl2"xy4" 
steel  plate 

M 

Mk23 

N 

81 -mm 
mortar 

P 

BDU-33 

Q 

Clutter 

R 

60-mm 

mortar 

S 

20-mm 

T 

30-mm 

U 

50-cal 

■ 
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Table  13  Description  of  target  buried  in  lanes  F~I, 


ID# 

Item 

Depth 

Azimuth 

Inclination 

XI 

60mm 

0.25 

0 

0 

X2 

60mm 

.5 

45 

30 

X3 

60mm 

.75 

90 

30 

X4 

60mm 

.75 

90 

45 

X5 

60mm 

1 

0 

90 

X6 

2.75” 

.25 

0 

30 

X7 

2.75” 

.5 

0 

45 

X8 

2.75” 

.5 

90 

60 

X9 

2.75” 

.75 

45 

75 

XIO 

81mm 

.25 

0 

30 

Xll 

81mm 

.5 

45 

45 

X12 

81mm 

.75 

90 

60 

X13 

BDU-33 

1 

45 

30 

X14 

1 20mm 

1 

45 

60 

NXl 

Steel  Bar  L24”/D3” 

.75 

0 

30 

NX2 

Steel  Bar 

.75 

90 

45 

NX3 

8”d  steel  plate 

.3 

Arbitrary 

0^ 

NX4 

8”d  steel  plate 

.5 

Arbitrary 

45^ 

NX5 

8”d  steel  plate 

.75 

Arbitrary 

90^ 

Scl 

Whole  crumpled  box  fin 

.25 

Arbitrary 

Arbitrary 

Sc2 

box  fin  fragment 

.5 

0 

0 

Sc3 

Projectile  nose 

.75 

0 

0l.3 

Sc4 

Cluster  of  clutter  items  25,17,  &  12 

.5^ 

0,120,240 

45^ 

Sc5 

clutter  Item  #22 

.75 

0 

30'’^ 

Sc6 

clutter  Item  #3 

.5 

0 

90^ 

Sc7 

clutter  Item  #24 

.75 

45 

30'’^ 

Sc8 

clutter  Item  #18 

.5 

90 

0 

Sc9 

clutter  Item  #23 

.75 

90 

03,5 

SclO 

clutter  Item  #11 

.5 

45 

OS 

O 

Sell 

clutter  Item  #16 

.5 

45 

45^ 

1 .  Curved  side  down 

2.  Depth  .5  m  to  upper  tip  of  eaeh  item,  tips  not  touching 

3.  No  twist  about  axis  after  inclination 

4.  45  deg  twist  after  inclination 

5.  Curved  side  up 
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